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ABSTRACT 


Two  concepts  for  the  augmentation  of  shock  tunnel  performance 
are  analyzed.  In  the  first  of  these,  termed  the  buffer  accelerator, 
a  magnetohydrodynamic  (MHD)  accelerator  is  used  to  accelerate  a 
buffer  gas,  which  in  turn  drives  a  shock  into  the  test  gas.  In  the 
second  approach,  termed  the  nozzle  accelerator,  an  MHD  accelera¬ 
tor  is  used  to  accelerate  the  test  gas  directly.  The  latter  approach 
is  shown  to  be  capable  of  providing  higher  velocities  and  longer  test 
times,  but  the  aerodynamic  quality  of  the  resulting  test  gas  suffers 
because  of  the  necessity  of  adding  seed  material  to  develop  the  con¬ 
ductivity  required  for  accelerator  operation.  Descriptions  are 
given  of  the  major  pieces  of  equipment  that  have  been  built  for  an 
experimental  program  to  investigate  the  two  augmentation  concepts. 
These  are  (1)  a  10^-joule  capacitor  bank  for  the  electric  field,  (2)  a 
10-weber/m^  coil  for  the  magnetic  field,  and  (3)  a  buffer  accelerator. 
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SECTION  I 
INTRODUCTION 


This  report  is  concerned  with  research  directed  toward  the  develop¬ 
ment  of  an  aerodynamic  test  facility  capable  of  providing  test-section 
conditions  corresponding  to  flight  in  the  atmosphere  at  velocities  from 
10,  000  to  36,  000  ft/sec,  at  altitudes  between  50,  000  and  250,  000  ft. 
Because  of  energy  considerations,  such  a  facility  must  operate  on  a 
pulsed  basis  rather  than  providing  a  continuous  flow.  A  companion 
report  (Ref.  1)  describes  related  work  on  a  shock  tunnel  system  which 
is  an  integral  part  of  the  effort. 

Impulse  aerodynamic  facilities  currently  operating  include  hotshot 
tunnels,  shock  tubes,  and  shock  tunnels.  The  limitations  on  the  per¬ 
formance  capabilities  of  hotshot  tunnels  are  fairly  well  established 
(Ref.  2),  Energy- density  limitations  of  the  arc  chamber  gas  limit  test- 
section  velocities  to  about  10,  000  ft/sec.  By  expanding  the  test  gas 
close  to  the  condensation  point,  these  tunnels  are  capable  of  providing 
flows  at  a  test- section  Mach  number  of  20,  at  unit  Reynolds  numbers 
corresponding  to  flight  at  about  100,  000  ft. 

Shock  tubes  are  capable  of  providing  stagnation  point  simulation 
(stagnation  pressure  and  stagnation  enthalpy  match  those  of  freeflight) 
for  the  conditions  of  interest  here.  Shock  tunnels  provide  Mach 
number-Reynolds  number  simulation  similar  to  that  provided  by  hotshots. 
The  limitations  on  shock  tube  performance  arise  from  a  number  of 
factors  (Refs.  1,  3,  and  4).  Arc  driven  tunnels  are  limited  by  the  energy 
loss  by  radiation  in  the  driver  (Ref.  5),  Another  limitation  is  the  practi¬ 
cal  maximum  pressure  that  can  be  contained. 

Based  on  the  results  from  hotshot  tunnels  and  shock  tunnels,  along 
with  a  limited  amount  of  free-flight  data,  dynamic  forces  and  heat- 
transfer  rat-cs  can  be  predicted  reasonably  well  from  the  data  from  these 
facilities.  Furthermore,  the  predicted  Mach  number  independence  of 
certain  coefficients  has  been  confirmed.  At  the  lower  densities  various 
second-order  effects  are  not  independent  of  Mach  number  but  depend 

upon  M/\/Re.  This  dependence  has  also  been  confirmed  (Ref.  6). 

At  the  higher  free- stream  velocities  of  interest  here,  it  is  expected 
that  Mach  number  independence  will  not  be  valid.  At  these  higher  energy 
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conditions,  the  finite  rates  of  the  various  reactions  taking  place  in  the 
gas  become  comparable  to  the  flow  transit  time.  To  provide  a  reason¬ 
able  simulation  of  these  effects,  it  is  necessary  to  match  free- stream 
velocities  (i.  e. ,  after-shock  enthalpy).  Matching  of  the  free -stream 
Mach  number  is  also  desirable  to  duplicate  the  flow'  field;  then  the 
free -stream  speed  of  sound  should  correspond  to  the  value  encountered 
in  freeflight.  To  duplicate  the  reaction  rates,  matching  of  the  pL  pro¬ 
duct,  where  p  is  the  free-stream  density  and  L  a  characteristic  length, 
should  also  be  imposed  (Ref.  7).  All  of  the  above  conditions  imply  that 
the  free-stream  Reynolds  number  is  also  matched.  However,  match¬ 
ing  of  the  pL  product  for  a  scale  model  implies  a  test- section  density 
higher  than  the  free-flight.  It  will  be  seen  that  obtaining  all  of  these 
conditions  simultaneously  is  extremely  difficult.  To  simplify  the  pre¬ 
sentation  of  performance  curves,  the  scale  effect  will  be  neglected, 
so  that  the  requirements  are  matching  of  velocity,  density,  and  Mach 
number. 

To  attain  the  conditions  as  given  above,  it  is  necessary  to  have  an 
airstream  of  high  velocity  and  high  density.  Because  of  the  limitations 
on  available  impulse  facilities  given  earlier,  it  is  apparent  that  a  new 
approach  is  necessary.  The  approach  chosen  here  is  to  augment  shock 
tube  performance  by  means  of  a  magnetohydrodynamic  (MHD)  accelera¬ 
tor  (Ref.  8),  in  which  directed  kinetic  energy  is  added  to  the  gas  by 
means  of  an  electrical  discharge  in  a  magnetic  field,  thereby  creating 
an  accelerating  body  force. 

The  MHD  accelerator  is  essentially  a  means  of  adding  energy  to  a 
moving  gas  stream.  Two  different  augmentation  approaches  are  con¬ 
sidered  here.  The  first  is  the  buffer  accelerator,  in  which  the  energy 
is  added  to  a  buffer  gas,  such  as  argon,  which  has  been  heated  by  a 
shock.  The  added  energy  causes  an  increase  in  shock  speed.  The 
shock  whose  speed  has  been  increased  is  then  driven  into  the  test  gas 
(air).  The  test  gas  is  then  expanded  to  the  test-section  conditions  by  a 
nozzle.  In  this  approach,  it  is  necessary  to  operate  in  a  nonreflected 
mode  (because  of  practical  pressure  limits),  with  the  result  that  test 
times  are  small  (fraction  of  a  millisecond).  However,  the  gas  in  which 
the  test  takes  place  is  theoretically  quite  clean. 

In  the  second  augmentation  approach,  termed  the  nozzle  accelera¬ 
tor,  the  MHD  accelerator  is  placed  in  the  nozzle  section  and  accelerates 
the  test  gas  directly.  This  is  a  much  more  efficient  process,  in  terms 
of  energy  addition,  than  the  buffer  accelerator.  However,  to  operate  the 
accelerator,  it  is  necessary  to  add  a  small  amount  of  easily  ionized 
seed  material  to  the  test  gas  to  develop  the  necessary  electrical  conduc¬ 
tivity  (Ref.  8).  The  effect  of  this  added  material  on  the  flow  about  the 
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test  object,  and  the  resulting  measurable  quantities  are  under  study.  It 
is  not  yet  clear  whether  this  seeding  will  present  a  major  problem. 

Following  a  section  on  general  considerations,  these  concepts  are 
discussed  in  some  detail  to  establish  the  performance  limits  attainable 
with  each.  The  remaining  sections  describe  equipment  that  is  being 
assembled  to  investigate  experimentally  the  two  concepts. 


SECTION  II 

SOME  GENERAL  CONSIDERATIONS 


There  are  some  considerations  associated  with  obtaining  a  high 
velocity,  high  density  flow  which  are  significant  regardless  of  the 
particular  approach  used;  they  are  discussed  in  this  section. 


2.1  ENERGY  REQUIREMENTS 

Figure  1  gives  the  values  of  the  flow  energy  in  a  usable  test  core 
of  7-ft  diameter.  Values  are  given  in  megajoules  per  millisecond  of 
run  time.  The  numbers  also  give  the  power  for  continuous  operation 
in  megawatts.  The  curve  of  100  megajoules  per  millisecond  corre¬ 
sponds  to  105  megawatts  on  a  continuous  basis.  This  is  approximately 
equal  to  the  total,  installed,  electric  generating  power  in  the  United 
States.  Therefore,  continuous  operation  in  this  regime  is  not  feasible. 
The  values  in  Fig.  1  are  highly  idealized  because  no  losses  have  been 
taken  into  account.  In  the  impulse-type  facilities,  the  input  energy  will 
be  higher  than  the  test- section  energy  by  a  factor  of  from  about  10  for 
highly  favorable  conditions  to  about  100  or  more  for  less  favorable 
conditions. 

The  attainable  run  time  will  probably  be  limited  to  about  1  msec. 

Run  time  limitations  result  both  from  the  limitations  of  the  gas  supply 
and  the  available  total  energy.  Existing  capacitor  energy  supplies  are 
limited  to  about  10  megajoules;  inductive  energy  supplies  can  supply  up 
to  100  megajoules,  but  probably  cannot  satisfy  the  1-msec  time  require¬ 
ment.  New  types  of  energy  sources  may  then  be  necessary.  Explosively 
driven  MHD  generators  may  fulfill  the  requirement. 


2.2  BOUNDARY  LAYERS  AND  HIGH  EXPANSION  RATIOS 

Boundary  layers  always  present  a  problem.  For  the  situation  con¬ 
sidered  here,  there  are  two  effects  of  interest,  both  of  which  cause  an 
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increase  in  the  total  energy  required.  The  first  of  these  is  the  boundary- 
layer  buildup  in  the  nozzle.  Figure  2  gives  estimates  of  the  geometric 
test- section  diameter  necessary  to  provide  a  specified  usable  core  diam¬ 
eter  at  the  indicated  altitude -velocity  point  {Ref.  9).  These  estimates 
are  based  upon  empirical  correlation  formulas  developed  for  lower  test- 
section  velocities.  To  have  a  usable  core  of  7  ft,  a  geometric  test- 
section  diameter  of  about  15  ft  is  required.  There  will  be  a  certain 
amount  of  flow  energy  in  the  boundary  layer  which  must  be  provided  for 
in  the  energy  source. 

The  second  effect  occurs  because  the  boundary  layer  places  a  lower 
limit  on  the  inside  diameter  of  shock  tubes  and  MHD  accelerators.  On 
the  other  hand,  large  area  expansion  ratios  are  required  to  reduce  the 
density  of  the  flow  to  the  desired  value.  The  minimum  shock  tube  diam¬ 
eter  and  the  necessary  expansion  ratio  may  combine  to  give  a  test- 
section  area  that  is  several  times  the  desired  value,  It  is  possible  to 
obtain  a  realistic  test-section  size  by  bleeding  off  part  of  the  test  gas, 
but  this  represents  a  loss  in  energy. 


2.3  RELAXATION  PROCESSES  AND  MATCHING  PROBLEMS 

In  a  nozzle  expansion,  the  static  temperature  falls  quite  rapidly. 

If  the  temperature  range  is  high  enough,  a  point  is  reached  at  which 
the  composition  of  the  gas  cannot  adjust  rapidly  enough  for  local  equilib¬ 
rium  to  be  attained.  The  composition  of  the  gas  in  the  test  section  will 
therefore  not  be  the  same  as  that  of  equilibrium  air  at  the  same  tem¬ 
perature  and  pressure. 

There  are  at  least  two  effects  of  this  mismatch  in  composition  that 
must  be  considered.  The  first  is  the  effect  of  the  nonequilibrium  com¬ 
position  on  the  flow  over  the  test  object,  and  the  influence  on  the 
resulting  measurable  quantities.  This  situation  is  not  yet  completely 
clear,  although  work  is  being  done  on  this  item  (see,  e.  g. ,  Ref.  10). 

The  second  factor  is  the  effect  of  this  altered  composition  on  the  nozzle 
flow,  and  these  effects  have  been  studied  extensively  (e.  g. ,  Refs.  11 
and  12).  There  is  first  of  all  a  reduction  in  the  velocity  at  the  nozzle 
exit,  since  a  portion  of  the  initial  thermal  energy  is  not  available  for 
conversion  to  kinetic  energy.  This  is,  for  the  application  here  con¬ 
sidered,  a  minor  effect.  More  important  is  the  alteration  in  the 
pressure- density  relation  caused  by  the  lower  temperatures.  At  a  point 
in  the  flow  at  which  the  temperature  is  a  given  value,  the  density  will  be 
higher  because  of  the  nonequilibrium  effects. 

Some  very  serious  questions  are  raised  by  the  fact  that  the  condi¬ 
tions  in  the  test  section  will  not  match  the  atmospheric  composition. 
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and  much  more  study  of  this  problem  is  needed  before  definite  answers 
are  possible.  For  an  equilibrium  expansion,  matching  of  density, 
velocity,  and  temperature  implies  matching  of  pressure  and  speed  of 
sound.  For  a  nonequilibrium  expansion,  this  is  not  true.  The  density 
and  velocity  may  be  matched  along  with  temperature  or  pressure  or 
speed  of  sound.  The  actual  choice  will  probably  depend  on  the  partic¬ 
ular  phenomenon  being  observed.  Here,  the  choice  is  made  to  match 
density,  velocity,  and  frozen  speed  of  sound  because  it  appears  that 
frozen  Mach  number  in  a  frozen  flow  most  closely  corresponds  to  the 
equilibrium  Mach  number  in  an  equilibrium  flow.  Further  study  may 
reveal  that,  for  some  conditions,  another  choice  should  be  made. 

In  recent  years,  much  has  been  learned  about  nozzle  processes, 
both  from  the  analytical  approach  of  Bray  (Ref.  11)  and  the  numerical 
studies  performed  at  Cornell  Aeronautical  Laboratory  (Ref.  12).  As 
yet,  sufficient  experimental  data  to  verify  these  results  have  not  been 
obtained.  For  the  conditions  of  interest  here,  the  primary  imbalance 
in  the  test  section  composition  is  in  the  atomic  oxygen  in  the  test  gas. 
Figure  3  gives  the  percentage  of  the  original  oxygen  that  is  in  atomic 
form  at  the  test  section  for  a  steady  isentropic  expansion.  It  should 
be  noted  that,  as  indicated  in  Fig.  3,  the  amount  of  the  atomic  oxygen 
frozen  in  expansion  is,  to  a  good  approximation,  independent  of  the 
reservoir  enthalpy  (or  test-section  velocity)  and  is  a  function  of  the 
entropy  alone.  For  the  condition  in  the  atmosphere,  the  entropy  is  a 
function  of  the  altitude  alone,  the  values  being  indicated  in  Fig.  3.  A 
value  of  the  parameter  f  of  0.  1  cm  corresponds  to  a  small  facility, 
whereas  a  value  of  10  cm  is  appropriate  for  the  large  facility  of  interest 
here. 

In  Fig.  4,  results  are  given  that  relate  altitude  to  entropy  for 
various  processes.  The  upper  curve  corresponds  to  an  equilibrium 
expansion,  *  the  middle  curve  to  an  expansion  with  chemical  freezing, 
with  matching  of  the  frozen  speed  of  sound,  and  the  lowest  curve  gives 
the  result  for  an  expansion  with  chemical  and  vibrational  freezing.  The 
degree  of  vibrational  freezing  is  currently  in  dispute,  but  the  differences 
are  not  a  major  factor  here.  The  middle  curve  is  arbitrarily  used  here. 

Thus  far,  no  way  has  been  found  to  overcome  the  effects  of  freezing. 
All  that  can  be  done  is  to  reduce  the  amount  of  mismatching  as  much  as 
possible  and  hope  that  corrections  to  the  experimental  results  can  be 
made  to  account  for  these  effects. 


Throughout  this  report,  atmospheric  properties  are  based  upon 
Ref.  13;  air  properties  are  from  Ref.  14  for  temperatures  above 
1500CK  and  from  Ref.  15  for  temperatures  below  1500°K. 
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2.4  HIGH  PRESSURES 

The  low  values  of  entropy  associated  with  the  low  altitudes,  as 
shown  in  Fig.  4,  indicate  that  very  high  equivalent  reservoir  pressures 
are  required.  The  pressure  level  will  therefore  have  to  be  placed  at 
the  highest  possible  value,  consistent  with  the  capabilities  in  pressure 
vessel  design. 


2.5  LOW  PRESSURES 

Although  very  high  pressures  are  required  at  one  point,  the  need 
exists  for  very  high  initial  vacuums  in  the  test  sections.  If  the  initial 
pressure  in  the  test  section  is  too  high,  a  secondary  shock  will  form 
during  the  nozzle  starting  process,  reducing  the  available  run  time 
(Ref.  16).  Figure  5  gives  the  initial  pressures  required  for  a  perfect 
start  of  the  nozzle  (no  secondary  shock).  Attainment  of  pressure  levels 
of  0.  1  and  0.  01  micron  of  mercury  will  require  a  great  deal  of  effort. 

It  should  be  noted  that  current  shock  tunnels  apparently  operate  at  initial 
pressures  of  10  or  100  times  that  required  for  a  perfect  start.  The 
reasons  for  this  are  beginning  to  be  evident  (see  Ref.  17),  but  as  yet, 
this  is  not  on  a  firm  basis. 

2.6  THE  DUPLICATION  LINE  AND  TESTING  AT  A  REDUCED  MACH  NUMBER 

In  describing  the  performance  capabilities  of  a  given  device,  use  is 
made  of  the  duplication  line  concept.  As  used  here,  duplication  means 
matching  of  density,  velocity,  and  frozen  Mach  number. 

For  a  given  physical  setup,  there  is  a  limiting  performance  line. 

The  limit  may  be  established  by  pressure  limitations,  energy  limitation, 
etc.  In  the  upper  part  of  Fig.  6,  such  a  limiting  duplication  line  is 
shown.  Duplication,  or  matching  of  density,  velocity,  and  Mach  num¬ 
ber,  is  possible  with  this  particular  setup  above  and  to  the  left  of  this 
line.  For  altitudes  above  that  indicated  by  A,  the  frozen  oxygen  content 
is  greater  than  that  at  A,  but  for  lower  altitudes  it  is  less. 

It  is  possible  to  test  at  lower  altitudes  than  those  above  the  limit 
duplication  line  by  relaxing  the  restriction  of  true  Mach  number  and  re¬ 
quiring  only  that  the  flow  be  hypersonic  (Mach  number  greater  than 
about  10).  Testing  at  a  reduced  Mach  number  with  matching  of  density 
and  velocity  involves  testing  at  a  higher  static  temperature  (i.  e.  ,  less 
expansion  of  the  nozzle  flow). 
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The  duplication  lines  such  as  shown  in  the  lower  part  of  Fig.  6  are 
generally  limit  lines  and  are  therefore  not  performance  lines  for  a 
given  physical  setup.  Determination  of  the  duplication  line  for  a  given 
physical  setup  is  a  more  difficult  task  and  is  not  discussed  herein. 

2.7  ENTROPY  PRODUCTION 

The  central  problem  of  providing  duplication  capability  in  the  de¬ 
sired  regime  can  be  related  to  entropy.  To  attain  the  desired  velocities, 
specific  amounts  of  energy  must  be  added  to  the  working  fluid.  To 
attain  the  altitude  duplication,  the  entropy  of  the  working  fluid  after  the 
energy  addition  must  not  exceed  a  certain  value.  Therefore,  the  energy 
must  be  added  very  efficiently,  in  terms  of  entropy  production. 

The  entropy  of  the  working  fluid  can  be  written  as  the  sum  of  two 
parts  -  the  initial  value  and  the  value  added  in  the  energy  addition  process 

S  =  Simiial  +  ASrnergv  addition 

Energy  addition  by  shock  heating  is  very  inefficient  in  terms  of  entropy 
production.  For  the  nonreflected  shock  tunnel,  it  is  necessary  to  work 
at  very  high  pressures  to  keep  the  initial  entropy  low.  For  the  MHD 
nozzle  accelerator,  the  entropy  production  can  be  controlled  to  a  certain 
extent,  but  the  requirement  for  a  low  value  of  the  initial  entropy  must  be 
compromised  with  the  conditions  of  temperature  and  pressure  necessary 
for  accelerator  operation. 


SECTION  III 

BUFFER  ACCELERATOR  PERFORMANCE 

3.1  DESCRIPTION  OF  OPERATION 

The  operation  of  a  shock  tunnel  with  a  buffer  accelerator  is  shown 
in  Fig.  7.  Diaphragm  D1  separates  the  high  pressure  driver  gas  from 
the  buffer  gas.  Upon  rupture  of  this  diaphragm,  a  shock  wave  SI  pro¬ 
pagates  into  the  buffer  gas.  The  plane  separating  the  two  gases  is 
termed  the  contact  surface,  *  Cl.  The  gas  behind  the  shock  SI  is  accel¬ 
erated  in  the  accelerator  so  that  the  shock  exiting  from  the  accelerator 
has  a  higher  velocity  than  it  did  before  entering  the  accelerator.  This 


Actually,  this  is  a  mixing  region  of  finite  thickness  rather  than  a 
plane. 
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higher  velocity  shock  then  ruptures  the  diaphragm  D2,  and  a  shock 
wave  S2  is  driven  into  the  test  gas.  The  gas  between  the  shock  S2  and 
the  contact  surface  C2  is  the  test  gas,  which  is  expanded  in  a  nozzle  to 
the  test- section  conditions.  (The  diaphragm  D3  is  inserted  to  allow  the 
nozzle  and  dump  tank  to  be  evacuated). 

This  concept  of  MHD  augmentation  in  the  buffer  of  a  shock  tunnel 
was  originally  proposed  by  Jan  van  der  Bliek,  ARO,  Inc.  ,  in  1963,  in 
an  unpublished  report.  Earlier,  work  had  been  done  by  Hogan  (Ref,  18) 
on  the  operation  of  an  accelerator  of  the  type  considered  herein.  His 
work,  which  was  continued  by  Leonard  (Ref.  19),  was  devoted  to  the 
understanding  of  the  operation  of  the  accelerator,  and  no  consideration 
was  given  to  driving  the  downstream  shock  into  a  test  gas.  A  device 
very  similar  to  the  one  considered  herein  is  being  studied  by  Leonard 
and  Rose  (Ref.  20)  in  which  the  accelerator  is  placed  at  the  diaphragm 
Dl,  The  driver  gas,  which  has  been  seeded,  is  accelerated,  causing  an 
increase  in  shock  speed  in  the  test  gas.  This  approach  involves  only 
two  sections  of  tube  rather  than  the  the  three  sections  that  are  required 
for  the  buffer  accelerator.  The  choice  between  these  two  devices  de¬ 
pends  upon  the  driver  capability.  Leonard  and  Rose  use  an  electric  arc 
driver  for  which  the  temperature  of  the  driver  gas  after  expansion  is 
sufficient  to  generate  the  necessary  conductivity.  At  VKF,  work  has 
been  devoted  to  drivers  with  indirect  heating,  for  which  the  temperature 
of  the  driver  gas  after  expansion  is  too  low  to  allow  for  operation  of  the 
accelerator.  Therefore,  the  shock-heated  buffer  gas  is  used  as  the 
working  gas  within  the  accelerator. 

3.2  ACCELERATOR  AS  A  SHOCK-ACCELERATING  DEVICE 

There  are  two  ways  of  studying  the  operation  of  an  accelerator  of 
this  type.  The  first  is  the  shock- acceleration  approach  of  Rosciszewski 
and  Oppenheim  (Ref,  21).  In  their  approach,  the  unsteady  x-t  charac¬ 
teristics  are  solved  to  obtain  shock  velocity  as  a  function  of  distance. 
They  are  able  to  obtain  an  approximate  solution  for  the  case  of  an  ideal 
gas  by  assuming  that  the  shock  is  a  U  +  a  characteristic.  This  approach 
has  not  been  followed  herein  because  extension  to  a  nonideal  gas  and 
removal  of  the  assumption  that  the  shock  follows  a  characteristic  would 
require  a  fairly  elaborate  computer  program. 

The  second  approach,  and  the  one  that  is  followed  herein,  is  the 
steady- state  approach.  It  is  assumed  that  the  accelerator  operation  is 
steady- state,  and  the  exit  conditions  may  be  obtained  by  applying  the 
steady- state  flow  equations.  Since  the  item  of  interest  for  the  buffer 
accelerator  is  the  downstream  shock  speed,  some  means  must  be  found 
for  relating  the  accelerator  exit  conditions  to  the  downstream  shock 
speed. 
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Figure  8  shows  the  situation  in  the  x-t  plane.  The  initial  shock 
is  accelerated  and  at  the  exit  condition  has  a  velocity  corresponding  to 
the  point  SA.  The  velocity  of  the  shock  at  this  point  is  that  obtained  by 
the  shock-acceleration  approach  of  Rosciszewski  and  Oppenheim.  A 
more  accurate  evaluation  requires  solving  the  characteristics  field  be¬ 
tween  the  shock  trajectory  and  the  indicated  U  -  a  wave.  After  the 
U  -  a  wave  has  traversed  the  accelerator,  steady- state  conditions  are 
established  in  the  accelerator.  The  shock  should  attenuate  after  it 
leaves  the  accelerator.  There  are  a  number  of  reasons  to  expect  this 
attenuation,  but  the  main  one  is  the  matter  of  adjusting  to  the  exit 
pressure  of  the  accelerator.  This  attenuation  of  the  shock  causes  the 
characteristics  to  coalesce  into  a  secondary  shock  which  moves  down¬ 
stream,  although  it  faces  upstream  (i.  e.  ,  the  downstream  pressure  is 
greater  than  the  upstream  pressure).  The  region  between  the  initial 
shock  and  the  secondary  shock  is  one  of  large  entropy  gradients  be¬ 
cause  adjacent  particles  have  passed  through  shocks  of  different 
strengths. 

An  alternate  interpretation  is  given  in  Fig.  9,  which  is  a  plot  of 
pressure  and  velocity.  The  shock  polar  gives  the  locus  of  conditions 
behind  a  shock  propagating  into  a  static  medium  with  a  pressure  Pj. 
Shock  speed  is  a  parameter  along  the  polar.  The  inlet  conditions  to 
the  accelerator  (denoted  by  i)  are  the  same  as  the  aftershock  condi¬ 
tions  (denoted  by  2)  at  the  inlet  shock  Mach  number.  The  shock  is 
accelerated,  and  the  conditions  at  the  exit  of  the  accelerator  at  the 
instant  that  the  shock  leaves  lie  at  the  point  SA  on  the  same  shock 
polar.  The  steady-state  velocity  at  the  accelerator  exit  (point  e)  is 
greater  than  the  inlet  velocity,  but  the  pressure  will  in  general  not  be 
too  different  from  the  inlet  pressure.  If  it  is  assumed  that  the  com¬ 
plicated  processes  are  adequately  described  by  simple  processes,  then 
the  simple  process  connecting  the  exit  conditions  e  with  the  shock  polar 
is  a  secondary  shock,  as  shown  in  the  figure  (Ref.  22,  p.  226).  The 
intersection  of  the  polar  for  the  secondary  shock  with  the  polar  of  the 
primary  shock  locates  the  points  b  and  c,  and  the  downstream  shock 
Mach  number  is  determined.  The  key  assumption  is,  of  course,  that 
the  processes  are  simple  ones.  This  should  be  a  fair  assumption,  since 
the  matching  is  done  on  the  basis  of  pressure  and  velocity,  which  are  not 
seriously  affected  by  the  entropy  gradients. 

This  secondary  shock  involves  a  loss  in  velocity,  and  therefore  its 
effect  should  be  minimized.  The  processes  are  quite  similar  to  the 
starting  processes  in  a  hypersonic  nozzle.  Based  upon  this  analogy, 

Dr.  W.  Warren*  has  suggested  placing  the  accelerator  next  to  the 


*General  Electric  Co.,  Space  Sciences  Laboratory,  Valley  Forge 
Space  Technology  Center,  King  of  Prussia,  Pennsylvania. 
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diaphragm  D2  and  lowering  the  initial  pressure  in  this  region  to  the 
point  that  no  shock  forms.  This  is  illustrated  in  Fig.  10.  The  practi¬ 
cal  disadvantage  of  this  approach  is  that  in  order  to  have  a  high  after¬ 
shock  pressure,  which  is  necessarj^  for  high  density  testing,  the 
accelerator  must  operate  at  high  pressure  at  the  exit  and  even  higher 
pressures  at  the  inlet. 

Thus,  far,  there  is  no  experimental  evidence  to  verify  the  existence 
of  the  secondary  shock.  The  available  experimental  data  (Refs.  18  and 
19)  indicate  a  large  deceleration  of  the  shock  at  the  accelerator  exit,  but 
the  deceleration  is  much  greater  than  can  be  attributed  to  the  formation 
of  the  secondary  shock.  This  large  deceleration  was  attributed  to  the 
rapid  reduction  in  the  magnetic  fields  at  the  accelerator  exit  (Refs.  18 
and  19). 

An  approximate  theory  can  be  developed.  Figure  11  illustrates  the 
nomenclature.  The  unsteady  processes  of  Fig.  8  have  been  omitted  for 
simplicity.  Ideal  gases  are  assumed.  Considering  the  downstream 
shock,  Mg^,  for  high  shock  strengths'  (Ms^  >  >  1), 


U„  =  Ms,/  a,/ 

(1) 

and 

fV  =  y\  ,  Ms,/  Pr 

(2) 

Across  the  secondary  shock. 

ue-u5/  j  2 

Pe  ‘ 

(3) 

a«  v  ^ 

/  +  Yi  ~  1 

1  Pt  Yx  +  l 

The  simplifying  assumption  is  now  introduced  that  (y^  - 

1)/(>1  +  1)  can 

be  neglected  in  comparison  with  P3  /Pe. 

Thus 

> 

y,  (y,  +  0 


(4) 


Across  the  interface  between  2'  and  3',  pressure  and  velocity  are 
matched. 

Pa'  =  P,'  (5) 

IV  =  IV  (6) 
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Combining  the  above  gives 


2>V 

v  +  1 


Pt 

p? 


Yx  (y,  +  i) 

yx*  (y,'  A  i) 


MSl. 


y  y  ,{v  +  1 )  /  P.. 

— 7-?r  m' mv  - 1 


0 


(7) 


This  is  a  quadratic  equation  which  gives  the  downstream  shock  speed, 
Mgj',  in  terms  of  the  exit  conditions,  Pe,  ae,  and  Me. 


Now  consider  the  case  for  which  there  is  no  diaphragm  at  the  exit 
of  the  accelerator,  so  that 


The  inlet  conditions  to  the  accelerator  are 


Pi  =  —r  MS;  p, 

Assume  that  the  accelerator  is  of  constant  area  so  that 

Pi  Uj  =  Ue 

and  assume  further  that  the  temperature  is  constant.  This  has  been 
found  to  be  a  reasonable  approximation  for  the  case  of  an  accelerator 
with  a  constant  area,  constant  electric  field,  and  constant  magnetic 
field.  With  these  assumptions,  the  following  result  is  obtained: 


Ms,.  _  i  jy  +  i 

Ms,  2  \j  y  —  l 

This  equation  gives  the  shock  speed  ratio  as  a  function  of  accelerator 
velocity  ratio.  Figure  12  gives  a  plot  of  Eq.  (8)  for  y  =  5/3.  Also 
shown  are  the  results  from  computer  solutions.  The  program  was 
written  by  Guy  Gilley,  ARO,  Inc. ,  and  it  solves  the  accelerator  equa¬ 
tions  in  the  form  given  in  Ref.  26  and  then  crosses  the  necessary 


ue 

+JO)  4- 

A  y~  J  +  A  > 

ly~  i 

Ui 

U, 

4  y  +  l  +  4  ' 

V  y  +  l 

Ue 

(8) 
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downstream  shocks.  Gas  properties  are  based  upon  a  singly  ionizing 
argon  model.  It  is  seen  that  Eq.  (8)  provides  a  reasonable  estimate. 


The  relationship  between  the  steady- state  accelerator  performance 
and  the  downstream  moving  shock  speed  given  by  Eq.  (8)  allows  the  in¬ 
crease  in  shock  speed  to  be  related  to  the  energy  added  in  the  accel¬ 
erator.  Assuming  that  most  of  the  energy  goes  into  the  gas  as  directed 
kinetic  energy,  *  then 


and  thus 


(1/2)  [ueJ 


E/Ai  -  (1/2)  p;  Ui* 


M  - 


(9) 


Furthermore,  the  strength  of  the  shock  in  the  test  gas  can  be  determined 
in  standard  ways.  Referring  to  Fig.  13,  either  an  unsteady  expansion 
fan  or  a  rear-facing  shock  occurs  between  region  b  and  region  (6). 

Given  the  above  results,  it  is  possible  to  relate  the  overall  perform¬ 
ance  to  the  energy  addition.  Figure  14  gives,  **  for  a  given  driver  (the 
600°K  externally  heated  driver  of  Ref.  1)  and  a  given  downstream  pres¬ 
sure,  Pfi,  the  energy  added  in  the  buffer  as  a  function  of  the  buffer  shock 
Mach  number,  Mg^  and  the  downstream  shock  Mach  number  Mgg.  It  is 
seen  that  for  a  given  downstream  shock  Mach  number,  there  is  an  opti¬ 
mum  value  for  the  buffer  shock  Mach  number.  The  optimum  shock  speeds 
are  high  enough  to  generate  the  conductivity  necessary  for  accelerator 
performance.  Using  the  minimum  energy  for  a  number  of  pressure  gives 
the  curves  shown  in  Fig,  15.  There  can,  therefore,  be  constructed  an 
equivalent  driver  line  as  the  basis  of  the  energy  added.  Such  a  set  of 
curves  is  given  in  Fig.  16,  For  reference  purposes,  the  line  for  an 
after-shock  pressure  (P7)  of  5000  atm  is  shown.  The  possibility  of  using 
helium  rather  than  argon  as  the  buffer  gas  has  also  been  investigated. 


& 

The  amount  of  energy  that  appears  as  thermal  energy  depends  upon 
the  specific  accelerator  design.  For  an  isothermal  accelerator  it  is 
zero,  but  for  the  constant  area  accelerator  with  constant  electric  field 
and  constant  magnetic  field  the  thermal  energy  addition  is  about  equal  to 
the  kinetic  energy  addition,  as  determined  by  computer  solutions. 

**These  results  were  obtained  from  a  computer  program  written  by 
Cecil  Taylor,  ARO,  Inc. 
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In  practice,  it  would  be  necessary  to  seed  the  helium  to  develop  the 
necessary  conductivity.  This  study  was  motivated  by  the  possibility  that 
the  energy  loss  by  radiation  from  the  shock- heated  argon  might  prove  to 
be  serious.*  The  results  indicate  that  helium  is  slightly  more  efficient 
than  argon  as  a  buffer  gas,  with  the  energy  added  being  about  2  5  percent 
less.  To  the  level  of  approximation  of  the  analysis,  this  amount  of 
energy  addition  is  not  particularly  significant. 


3.3  DETERMINATION  OF  TEST-SECTION  CONDITIONS 


Now  the  equivalent  shock  tube  performance  given  in  Fig.  16  will  be 
related  to  the  duplication  capabilities.  It  will  not  be  possible  to  reflect 
the  shock  wave  and  expand  through  a  converging- diverging  nozzle  because 
the  pressures  would  be  extremely  high  (20,  000  atm  or  greater).  There¬ 
fore,  the  expansion  of  the  gas  in  region  (?)  (Fig.  13),  begins  from  the 
shock  tube  area,  as  shown  in  Fig.  7. 


The  test- section  velocity  will  be  close  to  the  limiting  velocity 


where  *17^  is  the  stagnation  enthalpy  of  the  air  behind  the  shock  (region  (?)). 

Figure  17  gives  the  ratio  of  the  stagnation  enthalpy  to  the  static  enthalpy 
for  calculations  (Ref.  24)  based  upon  the  most  recent  air  data  (Ref.  14). 
For  simplicity,  this  ratio  will  be  taken  as  constant  at  1. 79.  The  varia¬ 
tions  about  this  number  are  minor.  Therefore, 

Vlim  ~\J( 2)  d-79)  ±-K 


From  Ref.  24,  Fig.  2d,  the  static  enthalpy  ratio  f^/hj  is  very  close  to 
the  ideal  gas  value,  which  for  y  -  1.4  and  large  shock  Mach  numbers 
is 


In  addition. 


so  that 


H,  =  3,234  x  lO*  ftVsec1 


(10) 


Vlim  =  1500  Ms, 


(ID 


*The  results  of  Ref.  23  indicate  that  at  the  high  pressures  of 
interest  herein  the  absorptivity  may  be  high  enough  that  there  will  be 
little  total  heat  loss,  in  which  case  this  will  not  be  a  problem. 
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with  V]^m  in  ft/sec.  To  obtain  the  entropy  of  the  gas  behind  the  shock, 
which  determines  the  duplication  altitude,  two  properties  are  needed. 
One  of  these  is  the  enthalpy,  which  is  given  by  Eq.  (10).  From  Ref.  24, 
Fig.  2a,  the  pressure  ratio  is  also  close  to  the  ideal  gas  value,  which 
is 


F S 
P. 


(12) 


With  the  pressure  and  enthalpy  determined,  the  entropy  can  be  obtained 
(Ref.  14)  and  the  altitude  determined  from  Fig.  4.  In  using  the  results 
of  Fig.  16,  the  shock  speeds  were  multiplied  by  0.  75  before  determining 
the  after- shock  conditions  to  make  an  approximate  allowance  for  shock 
attenuation. 


The  duplication  results  are  given  in  Fig.  18  for  the  equivalent  shock 
tube  performance  given  in  Fig.  16,  for  the  case  of  an  attenuation  factor 
of  0.  75.  It  is  seen  that  there  can  be  significant  improvement  in  the 
duplication  capabilities  by  the  addition  of  energy  in  the  buffer.  It  should 
be  pointed  out  that  very  large  nozzle  expansion  ratios  will  be  required. 


3.4  TEST  TIME  CONSIDERATIONS 

The  useful  test  times  with  this  device  will  be  quite  small.  Refer¬ 
ring  to  Fig.  19,  the  shock-tube  test  time  is 


lS.T.  -  U7 


Using 


gives 


u,  -  (i  - 


S'T-  ‘  MS  a 

Ss  ■  P* 

From  Fig.  2b  of  Ref.  24,  p^fpQ  is  about  8  for  the  conditions  of  interest 
here.  Therefore, 

l _ 


This  is  divided  by  two  to  account  approximately  for  the  test  time  loss 
attributable  to  interface  acceleration  (Ref.  25).  Thus, 

tS.T .  =  _14Ms  a 


There  is  a  loss  in  test  time  caused  by  the  convergence  of  a  U  -  a 
wave  from  the  throat  as  initiated  by  the  initial  shock,  and  a  U  +  a  wave 
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from  the  throat  as  initiated  by  the  contact  surface  (Ref.  16).  For  high 
test-section  flow  Mach  number  in  a  conical  nozzle. 


At 


L 


3.33 

MJ 


£t  ,s. 


The  flow  Mach  number  can  be  written  as 


M  = 

ftT.S. 

Using  Eq.  (11)  gives 


M 


1500 


a 


T.S. 


in  which  the  speed  of  sound  in  the  test  section,  arp  g  ,  is  in  ft /sec. 
Therefore, 


AtL 


3-33  L\ 
2^  V  aT.S. 


The  net  test  time  at  the  test  section  is  then 


lT.s.  =  ts.T.  “ 


L _ 3.33  lN 

14MSsas  2.25  MS,2  aT>s> 


For  simplicity,  let 

a,  =  a  t  $  =  1000  ft/sec  =  1  ft/mscc 


so  that,  with  the  test  time  in  msec  and  the  lengths  in  ft, 

^  _ _ L  _  3-33  Ljf 

T-s-  “  14  MS#  '  "  2_. 2 5  MSi* 


(15) 


The  shock  speed  for  which  the  test  time  is  maximum  for  given 
tube  length  and  nozzle  length  is  then 

Ms.  *  A*-  (16) 

and  at  this  condition 

Ms,  1 T  .s.  =  (17) 

If  the  tube  length  downstream  of  the  accelerator  is  100  ft,  and  the 
design  condition  such  that  M5g  =  20,  giving  a  test- section  Mach  num¬ 
ber  of  30,  then  the  test  time  is  180  ^sec  and  the  nozzle  length  48  ft 
At  lower  shock  Mach  numbers,  more  test  time  can  be  obtained  by 
proper  design,  but  the  result  is  that  test  times  will  generally  not  ex¬ 
ceed  200  fisec  (0.  2  msec).  It  should  be  again  noted  here  that  some 
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allowances  have  been  made  for  known  factors  affecting  test  time,  e.  g. , 
interface  acceleration  and  attenuation;  thus  the  above  test  time  esti¬ 
mates  are  as  realistic  as  possible  based  on  the  present  state-of-the-art. 


3.5  DETAILS  OF  ACCELERATOR  OPERATION 


In  the  preceding  analysis,  the  accelerator  was  treated  as  an  energy 
addition  device,  but  no  consideration  was  given  to  the  details  of  the 
accelerator.  Now  some  basic  considerations  will  be  applied  to  determine 
the  size. 


The  one- dimensional,  steady-flow  equations  describing  the  acceler¬ 
ator  operation  are  (Ref.  26) 

Conservation  of  mass 

pUA  =  m  =  constant  (18) 


Conservation  of  momentum 


pU 


du  ,  dP 
d  x  dx 


exB (E  -  BU) 


Conservation  of  energy 


pU 


crE(E  -  BU) 


(19) 

(20) 


(The  nomenclature  is  illustrated  in  Fig.  20.  )  The  system  of  equations 
is  completed  by  the  equation  of  state  and  the  specification  of  three 
quantities.  For  the  buffer  accelerator,  it  is  reasonable  to  assume  that 
E,  B,  and  the  cross-sectional  area.  A,  are  constants.*  A  further 
reasonable  approximation  is  to  neglect  the  pressure  gradient  term  in 
Eq.  (19)  and  assume  that  the  conductivity  is  constant.  This  equation 
may  then  be  integrated  to  give 


The  asymptotic  final  value  if  E/B,  and  this  velocity  is  achieved  in  a 
length  for  which  the  interaction  parameter 


is  about  2  or  3. 


S 


oB1  La 
pU 


(22) 


*It  may  be  necessary  to  have  some  variation  in  E  to  overcome  the 
maximum  velocity  increase  limitation  pointed  out  in  Ref.  2  6. 
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To  illustrate  the  relative  magnitudes  involved,  the  condition  of 
Msg  =  20,  Pg  =  10  atm  is  chosen.  The  following  values  were  obtained 

as  part  of  the  computer  solutions  upon  which  Figs.  14,  15,  and  16  are 
based: 

Pt  =  1.7  atm 
IWUi  =  8.5 
Ms,  =  11 

From  these,  the  desired  value  of  E/B  is 

=  24,800  m/sec 

B 

The  estimated  conductivity  is  5  x  10*3  mho/m,  and  assuming  an  inter¬ 
action  parameter  of  3  gives 

B2  La  =  23 

in  which  B  is  in  weber/m^,  and  the  length  La  in  meters.  Therefore, 
an  accelerator  1  m  long  with  a  B-field  of  5  weber /m2  is  adequate.  The 

E- field  is  then 


E  -  124,000  v/m  =  1240  v/cm 


For  an  electrode  separation  distance  of  20  cm,  the  electrode  voltage 
is  24,  800  v.  From  Fig.  15,  the  energy  flow  rate  is  5  x  108  joule/ 
sec-cm^  so  that  for  a  400-^sec  run*  and  a  cross-sectional  area  of 
400  cm2,  the  ideal  total  energy  required 


Energy  =  (5  x  10s)  (400  x  10“‘)(400)  =  8  x  107  joules 


The  average  total  current  is  then 

T  =  _ 8  *--107  .... 

(400  x  10“ 4 ) ( 24,800 ) 


8  x  106  amp 


This  is  quite  a  high  value.  ♦*  The  current  density  is  concentrated  at  the 


♦The  accelerator  run  time  should  be  longer  than  the  shock  tube  test 
time  to  allow  for  the  fact  that  disturbances  from  the  accelerator  travel 
at  the  speed  of  sound  relative  to  the  flow.  The  amount  of  this  has  not 
been  estimated. 

♦♦The  high  value  of  the  total  current  suggests  the  possibility  of  using 
this  current  to  generate  the  magnetic  field  by  means  of  current  loops  in 
series  with  the  electrodes. 
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leading  edge,  and  the  initial  value  is  found  to  be 

Jmax  =  5.45  x  10*  amp/cms 
This  is  also  a  high  value. 


3.6  SUMMARY  OF  BUFFER  ACCELERATOR 

These  estimates  have  shown  that  the  buffer  accelerator  is  capable 
of  providing  the  capability  for  testing  in  a  region  which  is  not  currently 
available.  In  the  next  section,  the  nozzle  accelerator  is  considered, 
and  it  will  be  seen  that  this  approach  is  also  theoretically  capable  of 
providing  coverage  of  a  somewhat  similar  region.  In  many  respects, 
the  nozzle  accelerator  is  considered  to  be  superior  to  the  buffer  accel¬ 
erator.  However,  the  buffer  accelerator  has  one  major  advantage: 
the  test  gas  is  clean.  For  the  nozzle  accelerator,  it  is  necessary  to 
add  seed  material  to  the  test  gas;  therefore  there  is  foreign  material 
in  the  test  gas.  The  choice  between  a  buffer  accelerator  and  a  nozzle 
accelerator  lies  wholly  in  the  seriousness  with  which  the  effect  of  the 
seed  addition  is  viewed. 

The  primary  disadvantage  of  the  buffer  accelerator  is  the  short 
test  time,  which  is  on  the  order  of  200  ^sec.  It  may  be  possible  to 
perform  meaningful  aerodynamic  experiments  within  this  time  interval. 

The  success  of  a  buffer  accelerator  depends  upon  a  number  of 
factors.  The  test  program  planned  for  the  buffer  accelerator  has  the 
following  purposes: 

1.  Determination  of  the  validity  of  the  mathematical  models  used, 

2.  Determination  of  the  overall  efficiency,  and 

3.  Investigation  of  the  ability  to  draw  large  current  densities  in 
a  diffuse  sheet. 


SECTION  IV 

NOZZLE  ACCELERATOR  PERFORMANCE 

4.1  DESCRIPTION  OF  OPERATION 

A  simplified  schematic  of  a  nozzle  accelerator  is  shown  in  Fig.  21. 
The  shock  tube  (which  need  not  have  a  buffer  section)  generates  a  slug 
of  high  pressure,  high  temperature  air.  *  The  shock  tube  may  be 


*An  alternate  approach  is  to  use  a  hot  shot  as  the  gas  source  (Ref.  28). 
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operated  in  a  nonreflected  mode,  but  in  general  the  requirements  of  run 
time  will  dictate  operation  in  the  reflected  mode,  using  either  the 
tailored-interface  or  the  equilibrium  interface-technique  (Ref.  27).  The 
test  gas  from  the  shock  tube  is  then  accelerated  by  the  MI1D  accelerator. 


4.2  DISCUSSION  OF  ADVANTAGES  AND  DISADVANTAGES 

The  primary  advantage  of  this  approach  is  that  the  process  is  in¬ 
herently  more  efficient  than  the  buffer  accelerator  because  the  energy 
is  added  directly  to  the  test  gas.  The  disadvantage  is  the  necessity  of 
adding  seed  material  to  the  test  gas,  since  at  the  conditions  for  which 
air  has  any  significant  electrical  conductivity  the  entropy  exceeds  that 
necessary  for  flight  duplication  at  the  altitudes  of  interest. 

The  intentional  introduction  of  a  foreign  material  into  the  test 
gas*  must  be  viewed  with  suspicion  especially  because  of  the  large 
effort  that  was  necessary  to  reduce  the  contamination  of  hotshot  tun¬ 
nels  to  the  point  where  the  data  were  useful  (Ref.  2).  However,  seed 
is  necessary  for  the  nozzle  accelerator  to  be  a  workable  device.  The 
effects  of  this  seed  material  upon  the  flow  about  the  body  are,  at 
present,  uncertain  although  the  problem  is  being  studied  (Ref.  30). 

Some  order-of-magnitude  estimates  given  in  Ref.  31  indicate  that  the 
seed  may  not  adversely  affect  the  flow,  particularly  at  velocities  above 
30,  000  ft/sec. 

4.3  DUPLICATION  CAPABILITIES 

This  particular  approach  has  received  extensive  study.  The  theory 
is  fairly  well  developed,  and  there  exists  a  general  optimization  theory 
(Ref.  31).  Other  analytical  studies  are  given  in  Refs.  8,  32,  33,  and 
34,  which  may  be  taken  as  representatives  from  a  large  body  of  avail¬ 
able  papers  and  reports.  Experimentally,  the  emphasis  has  been  upon 
accelerators  for  continuous  operation  (see  Ref.  35),  but  the  possibility 
of  use  with  pulsed  supplies  is  receiving  increased  attention  (Refs.  28, 
36,  37,  and  38). 

Because  a  considerable  body  of  literature  exists  for  the  nozzle 
accelerator  approach,  no  attempt  will  be  made  to  present  the  analytical 
details.  Instead,  results  from  previous  work  will  be  cited. 


A  means  of  seed  injection  is  given  in  Ref.  29. 
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The  particular  results  cited  here  are  those  given  in  Ref.  32.  The 
analysis  is  based  upon  a  set  of  operating  conditions  chosen  as  follows: 

1.  The  applied  magnetic  field  (B)  is  constant, 

2.  The  static  temperature  is  constant,  and 

3.  The  loading  factor  (E/BU)  is  constant. 

The  inlet  conditions  are  chosen  to  maximize  the  test- section  velocity 
for  a  given  entropy.  Although  the  particular  set  of  operating  condi¬ 
tions  was  chosen  on  the  basis  of  analytical  simplicity,  the  results  from 
this  approach  agree  well  with  those  of  a  more  exact  theory,  in  which 
the  optimum  operating  conditions  are  determined  (Ref.  32).  In  addition, 
studies  have  been  made  of  accelerators  with  some  practical  simplifica¬ 
tions  (primarily  in  the  area  distribution). 

These  studies  show  that  the  duplication  limits  are  determined  by 
three  basic  parameters: 

1.  seed  concentration  (i.  e. ,  conductivity); 

2.  the  parameter  B2La,  where  B  is  the  externally  applied  mag¬ 
netic  field  strength  and  La  the  accelerator  length;  and 

3.  the  stagnation  pressure  of  the  reservoir  supplying  the 
accelerator. 

A  high  value  for  each  of  these  is  called  for.  The  first  two  of  these 
parameters  should  actually  be  considered  together  since  the  governing 
combination  is  really  the  interaction  parameter  (Eq.  22))  which  con¬ 
tains  the  grouping  <rB2La.  The  conductivity  varies  roughly  as  the  square 
root  of  the  seeding  fraction;  thus,  calculations  based  upon  one  value  of 
the  seeding  fraction  can  be  used  for  other  seeding  fractions  by  making  the 
proper  adjustment  in  the  parameter  B2La. 

The  requirement  for  a  high  seeding  fraction  must  be  compromised 
with  the  difficulty  of  data  interpretation  if  the  test  gas  contains  a  large 
amount  of  seed.  Since,  at  present,  the  degree  of  difficulty  caused  by 
seed  is  uncertain,  the  allowable  seed  fraction  must  be  set  arbitrarily. 

The  studies  that  have  been  conducted  to  date  have  been  based  upon  either 
1/4-  or  1 /2-percent  potassium  seed.  The  actual  contaminant  concentra¬ 
tion  will  probably  be  larger,  since  the  potassium  would  in  practice  be 
introduced  in  the  form  of  a  compound. 

The  upper  limit  on  the  parameter  B2La  may  be  set  by  a  number  of 
factors.  At  very  high  B-field  levels,  the  internal  magnet  stresses  may 
be  large.  At  a  magnetic  field  strength  of  25  weber/m2  this  stress  is 
36,  000  psi.  By  proper  choice  of  material,  this  may  be  contained.  Joule 
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heating  is  a  problem  for  continuous  operation,  but  it  is  probably  not  a 
factor  here  since  the  coil  may  be  operated  on  a  pulsed  basis.  Further 
discussion  is  given  in  Section  VII.  The  length  of  the  accelerator  must 
not  be  such  that  boundary -layer  buildup  within  the  accelerator  is  a 
problem.  Other  practical  problems  may  also  limit  the  length.  The 
energy  in  the  B  field,  and  therefore  the  energy  required  to  generate  the 
field  on  a  pulsed  basis,  is  proportional  to  the  product  of  B2  and  the 
volume  over  which  the  field  acts;  therefore  the  energy  varies  directly 
as  B2La.  The  actual  energy  necessary  can  be  determined  only  by  a 
detailed  design,  taking  into  account  the  various  packaging  problems. 

The  results  of  the  calculations  are  given  in  Figs.  22  and  23.  Re¬ 
sults  are  given  for  three  different  values  of  the  stagnation  pressure, 

100,  1000,  and  5000  atm,  representing  three  different  levels  of  shock- 
tube  performance.*  The  influence  of  the  accelerator  performance  is 
shown  by  the  three  values- of  B2La  of  100,  500,  and  2000  weber2/m3. 

The  first  of  these  values  corresponds  to  a  magnetic  field  of  5  weber /m2 
with  an  accelerator  length  of  4  m  {or  other  possible  combinations).  The 
higher  values  correspond  to  accelerator  lengths  of,  say,  5  m,  with  re¬ 
sulting  magnetic  fields  of  10  weber/m2  and  20  weber/m2,  respectively. 

The  figures  illustrate  first  of  all  that  this  type  of  accelerator  is 
limited  to  duplication  of  altitudes  above  150,  000  to  200,  000  ft.  In 
Ref.  8,  the  limits  are  given  as  from  100,  000  to  150,  000  ft.  There  are 
two  reasons  for  these  differences.  The  first  is  that  here  the  more  re¬ 
strictive  requirement  of  true  test-section  Mach  number  is  imposed, 
rather  than  the  requirement  of  a  hypersonic  Mach  number  used  in  Ref.  8. 
The  second  reason  is  that  the  relationship  between  altitude  and  entropy 
used  in  Ref.  8  is  based  upon  an  expansion  with  frozen  chemistry  and 
vibration,  whereas  here  the  relationship  is  based  upon  an  expansion  with 
frozen  chemistry  and  vibrational  equilibrium. 

Figures  22a,  b,  and  c  show  the  effect  of  higher  accelerator  per¬ 
formance  for  a  given  supply  stagnation  pressure.  The  results  show 
quite  strongly  the  necessity  of  high  values  of  B2La  if  high  test-section 
velocities  are  to  be  obtained  at  altitudes  below  250,  000  ft.  For  a  supply 
pressure  of  100  atm,  the  curves  are  terminated  on  the  low-velocity  end 
at  the  velocity  for  which  the  inlet  Mach  number  is  1.5.  The  low- velocity 
termination  point  for  the  case  of  a  5000-atm  supply  is  the  limit  of  the 
tables  of  Ref.  19. 


♦Shock  tube  performance  depends  upon  stagnation  temperature  as 
well  as  stagnation  pressure,  but  for  these  calculations  the  temperature 
is  implicitly  determined  by  the  requirement  of  adequate  electrical 
conductivity. 
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Figure  23  illustrates  the  effect  of  the  supply  for  a  given  accelerator 
(B^La  =  500  weber^/m^).  From  this  it  is  seen  that  the  influence  of  the 
stagnation  pressure  is  not  as  strong  as  that  of  the  accelerator  perform¬ 
ance,  which  is  perhaps  not  too  surprising. 

In  comparing  the  capabilities  of  the  nozzle  accelerator  as  given  in 
Fig.  22  with  those  of  the  buffer  accelerator  as  given  in  Fig.  18,  it  is 
seen  that  the  nozzle  accelerator  provides  a  capability  at  the  high 
velocities  that  the  buffer  accelerator  cannot  provide.  The  buffer  accel¬ 
erator,  on  the  other  hand,  provides  coverage  of  the  high  density  (low 
altitude),  moderate  velocity  area,  which  the  nozzle  accelerator  does 
not  cover. 

4.4  SPECIFIC  CASE 

To  illustrate  the  magnitudes  of  the  parameters  involved,  a  specific 
case  will  be  assumed.  Take  from  Fig.  22b  the  case  for  a  test-section 
velocity  of  30,  000  ft/sec  at  a  duplication  altitude  of  200,  000  ft.  The 
parameter  B^La  is  then  2000  weber^/m^,  which  could  be  provided  by 
a  magnetic  field  strength  of  20  weber/m^,  and  an  accelerator  length  of 
5  m.  *  A  shorter  accelerator  with  larger  magnetic  field  strengths  would 
be  desirable.  The  reservoir  pressure  is  1000  atm. 

From  Ref.  8,  the  optimum  accelerator  inlet  conditions  are 

T;  =  3000° K 
P.  =  7.67  atm 
S/R;  =  31.4 

Because  of  the  uncertainty  in  the  conductivity  estimates,  the  inlet  condi¬ 
tions  are  taken  at  a  higher  temperature.**  Therefore,  take 

Tj  =  3500°K  S/Rj  =  31.4 

P;  =  21.0  atm  Pi  =  0.129  Ib/fP 


*Smaller  values  are  obtained  if  the  seeding  fraction  is  increased 
above  1/4  percent. 

**Recent  unpublished  studies  have  shown  that  in  order  to  keep  the 
magnitude  of  the  Hall  potential  reasonable,  even  higher  inlet  tempera¬ 
tures  are  required. 
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The  reservoir  conditions  are  determined  by  the  stagnation  pressure 
and  the  inlet  entropy 

P0  =  1000  atm 
S/Ro  =  31.4 

from  which  T0  =  6100°K 

The  shock  tube  conditions  required  to  obtain  these  conditions  in  the  re¬ 
flected  region  are 

Ms  =  8.8 

Pj  =  1  atm 

Therefore 

Uj  =  10,450  ft/sec 
M.  =  2.9 

The  mass  flow  rate  through  a  usable  core  of  7 -ft  diameter  is 
m  *  p  AU 

=  (1.29  X  10" S)  (38.5)  (30,000) 

=  15  lb /sec 

From  Fig.  2,  the  geometric  diameter  would  have  to  be  about  15  ft.  The 
mass  flow  entrapped  in  the  boundary  layer  in  the  test  section  is  not 
definite,  but  the  total  mass  flow  rate  will  be  less  than  four  times  that 
through  the  core.  Therefore,  as  an  approximation,  take  the  mass  flow 
rate  through  the  accelerator  as 

m  =  50  Ib/sec 

This  assumes  that  the  boundary- layer  buildup  occurs  primarily  down¬ 
stream  of  the  accelerator.  The  inlet  area  is  then 

A.  =  m/p;Ui 

=  50/(0.129)  (10,450) 

=  0.0371  ft2 
=  5.34  in.2 

which  requires  an  inlet  diameter  of  2.  31  in. 


23 


A  £  DC-  T  R-66-25 


The  accelerator  exit  conditions  are  found  to  be 

Te  =  3500° K 

Pe  =  P./e  7.73  atm 

pe  =  p./e  =  0.0475  lb/fts 

<S/R)e  =  32.4 

U£  =  29,800  ft/sec 

Me  =  8.2 

Ae  =  5.1  in.J 

An  exit  diameter  of  2.  26  in.  is  required.  Continuity  considerations 
therefore  call  for  a  slight  decrease  in  flow  area  through  the  accelerator, 
which  is  undesirable  for  other  reasons.  In  practice,  a  sufficient  in¬ 
crease  in  area  must  be  provided  to  eliminate  the  possibility  of  choking 
caused  by  boundary- layer  formation.  The  amount  necessary  is  not 
known  absolutely,  but  it  is  assumed  that  a  small  additional  divergence 
angle  of  1  or  2  deg  on  the  walls  of  the  accelerator  will  be  sufficient. 

By  use  of  a  real  air  expansion,  the  nozzle  throat  area  is  found  to  be 
0.  00805  in.  2,  calling  for  a  diameter  of  0.  2  8  in. 

Turning  now  to  the  power  supply  requirements  for  the  accelerator, 
the  energy  added  in  the  accelerator  is  ideally 

E  -  if.  (Ue5  -  Uj’1/2 

=  8  x  10'  joule/sec 

=  0.8  megajoule/ir.sec 

The  power  supply  requirement  then  depends  upon  the  run  time  and  the 
overall  efficiency.  A  run  time  of  2  or  3  msec  should  be  possible  at  this 
shock  tube  condition.  The  overall  efficiency  is  unknown,  but  assuming 
provisionally  a  value  of  10  percent  gives  as  the  energy  requirement  or 
24  megajoules  for  3  msec. 

The  maximum  voltage  required  may  be  estimated  by  using  the  fact 
that  at  the  exit  the  loading  factor 

£  =  E/BU 

is  close  to  unity.  Therefore 

E  =  (20)  (29,800  x  0.3048) 

=  182,000  v/m 
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The  po  ential  difference  across  the  core  is 

V  =  (182,000)  (2.26/39.37) 

=  10,450  v 

To  this  must  be  added  the  voltage  drop  in  the  boundary  layers,  which 
is  at  present  unknown.  In  any  event,  the  E-field  power  supply  must 
have  a  voltage  capability  of  from  10  to  2  0  kv. 

The  accelerator  will  be  built  of  segmented  electrodes,  with  indi¬ 
vidual  floating  power  supplies  for  each  pair,  to  neutralize  the  Hall 
voltage  (Ref.  8).  The  magnitude  of  this  Hall  potential  is  estimated  at 
from  10  to  40  kv.  There  are  a  number  of  problems  to  be  studied  in 
regard  to  the  number  of  electrodes  required  and  the  insulation  thick¬ 
ness  between  adjacent  electrodes  that  is  necessary  to  prevent  shorting 
caused  by  the  Hall  potential. 

Estimates  of  the  energy  requirement  for  the  magnetic  field  power 
supply  are  more  difficult  to  obtain  since  the  total  volume  over  which 
the  B-field  is  generated  must  be  known.  This  total  volume  must  take 
into  account  the  increase  in  accelerator  size  attributable  to  boundary- 
layer  growth  as  well  as  the  various  packaging  problems.  At  a  B-field 
level  of  20  weber/m^,  the  energy  stored  in  the  field  is  about  8x  10**  joules 
per  cubic  meter  of  unit  inside  coil  volume.  The  internal  volume  of  the 
accelerator,  not  including  boundary-layer  corrections,  is  0.  0168  m^. 

On  this  basis,  the  ideal  energy  requirement  is  1.34  x  10?  joules.  There¬ 
fore,  the  actual  energy  requirement  to  provide  the  B  field  on  a  pulsed 
basis  will  be  perhaps  100  megajoules.  This  is  indeed  a  large  value,  and 
study  must  be  devoted  to  this  particular  aspect  of  the  problem  to  keep  it 
within  bounds.  It  is  interesting  to  note  that  the  energy  requirement  to 
provide  the  B-field  is  greater  than  that  required  to  provide  the  E-field. 

The  purpose  in  presenting  this  specific  case  is  not  to  provide  an 
actual  design  but  to  point  out  the  problem  areas  and  areas  of  uncertainty. 
Design  of  an  actual  facility  requires  answ'ers  to  some  of  these  problems 
as  well  as  more  complete  parameter  studies. 


4.5  SUMMARY  OF  NOZZLE  ACCELERATOR 

The  nozzle  accelerator  has  the  theoretical  capability  of  providing  a 
test  capability  approaching  escape  velocity  for  altitudes  above  150,  000  to 
200,  000  ft. 
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The  primary  areas  that  need  study  are 

1.  seeding,  both  in  regard  to  technique  and  effect  on  flow  about 
test  object; 

2.  boundary  layers,  both  in  the  accelerator  channel  and  the  nozzle; 

3.  determination  of  overall  efficiency; 

4.  working  with  high  voltages  and  high  energies  on  a  pulsed  basis; 
and 

5.  electrode  segmentation. 


SECTION  V 

BASIC  APPROACH  FOR  EXPERIMENTAL  STUDIES 

The  two  previous  sections  have  pointed  out  that  of  the  two  possible 
approaches  to  MHD  augmentation  of  a  shock  tube,  the  buffer  accelerator 
theoretically  provides  a  clean  test  gas  but  short  test  times,  whereas  the 
nozzle  accelerator  requires  the  addition  of  a  seed  material  to  the  test 
gas.  The  choice  between  the  two  depends  upon  the  relative  importance 
of  a  clean  test  gas  and  test  time  for  the  class  of  problems  to  be  studied. 
For  aerodynamic  tests  the  clean  test  gas  will  probably  be  the  governing 
consideration,  whereas  for  propulsion  tests  the  test  time  may  be  the  key 
item. 

Because  of  the  lack  of  data  for  either  approach,  it  was  decided  to 
perform  experimental  studies  for  both.  Following  these  preliminary 
studies,  one  of  the  two  approaches  will  be  selected  as  the  most  prom¬ 
ising,  and  further  work  will  be  performed. 

The  guiding  consideration  in  the  present  research  program  is  to 
develop  a  test  setup  that,  with  minor  modification,  can  be  used  for 
investigating  both  the  buffer  accelerator  and  nozzle  accelerator  con¬ 
cepts.  For  either  of  these,  there  are  four  basic  items  needed:  (1)  a 
shock  tube,  (2)  an  electric  field  supply,  (3)  a  magnetic  field  power 
supply  and  a  magnet,  and  (4)  an  accelerator.  Of  these  four  items,  the 
first  three  are  common;  the  major  difference  is  in  the  accelerators. 

The  buffer  accelerator  can  consist  of  a  single  pair  of  electrodes, 
whereas  the  nozzle  accelerator  will  have  a  number  of  segmented  elect¬ 
rodes.  This  difference  leads  to  certain  differences  in  the  power  supply 
requirements,  and  of  necessity  there  will  be  modifications  to  the  tube 
to  accommodate  the  accelerators.  Fortunately,  these  differences  are 
not  major. 

In  discussing  the  components  in  the  remainder  of  the  report,  the 
design  for  the  buffer  accelerator  is  given.  The  buffer  accelerator  was 
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chosen  for  initial  investigation,  because  an  AEDC-funded  program  on 
the  nozzle  accelerator  was  underway  elsewhere  (Ref.  36). 

Considering  first  the  shock  tube,  the  existing  Tunnel  J  shock  tube 
was  selected  because  of  its  availability  and  performance.  Planned 
development  of  drivers  for  this  tube  is  adequate  for  the  purposes  here. 
Reference  1  gives  in  detail  the  capabilities  and  planned  development  of 
this  tube. 

The  electric  field  power  supply  must  be  sufficiently  flexible  to 
operate  as  a  single  bank  for  the  buffer  accelerator  and  as  a  number  of 
independent  banks  for  the  nozzle  accelerator.  Capacitors  fulfill  this 
requirement  quite  well;  therefore  the  existing  106- joule  Hotshot  1  capa¬ 
citor  bank  (Ref.  2)  was  relocated.  This  particular  bank  consists  of 
4000-v  cans,  and  stacking  is  necessary  to  develop  the  required  voltages. 

A  description  of  this  bank  is  given  in  Section  VI. 

The  decision  was  made  to  use  the  existing  Tunnel  F  generators  as 
the  power  supply  for  the  magnetic  field  and  to  procure  a  special  coil. 

Since  the  generators  are  of  relatively  low  voltage,  a  large  cross- 
sectional  area  in  the  coil  is  required.  By  judicious  design,  a  potential 
magnetic  field  strength  of  10  weber/m^  was  obtained,  which  is  sufficiently 
high  to  perform  meaningful  experiments.  The  possibility  also  exists  for 
increasing  the  attainable  field  strength.  The  details  of  this  component 
are  given  in  Section  VII. 

The  final  component  is  the  accelerator.  The  design  for  the  buffer 
accelerator  is  given  in  Section  VIII.  The  design  of  the  nozzle  accelerator 
is  not  complete. 


SECTION  VI 

ELECTRIC  FIELD  SYSTEM 


6.1  INTRODUCTION 

Because  of  considerations  presented  in  previous  sections  of  this 
report  and  other  factors,  the  general  design  criteria  for  the  electric 
field  system  were  to  deliver  10®  joules  of  electrical  energy  into  a  con¬ 
ditioned  gas  behind  the  shock  wave  within  a  specific  time  duration,  with 
a  given  pulse  shape,  and  at  specific  electric  field  and  current  density 
levels.  The  conductivity  of  the  accelerator  inlet  gas  could  be  varied 
one  and  one-half  orders  in  magnitude,  depending  on  the  particular  shock 
tube  performance  regime. 
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The  specific  electric  field  system  design  criteria  were  determined 
to  be  as  follows: 

(a)  Total  Energy:  10®  joules 

(b)  Gas  Conductivity:  10®  -  5  x  104  mho/m 

( c )  Electric  Field  Levels:  0-10**  v/m  (0-4000  v) 

(d)  Current  Density  Levels:  0-4  x  10®  amp/m 

(B  =  0)  (0  to  1.  3  x  10®  amp) 

(e)  Pulse  Time  Duration:  400-2000  ftsec 

(f)  Pulse  Shape:  Square 

(g)  Electrode  Configuration:  Internal  surfaces  to  match  interior 
of  40-mm  tube  and  of  variable  length  from  1/2  to  12  in. 

(h)  Pulse  Synchronization:  Microsecond  resolution  in  firing 

pulse  relative  to  shock  wave 

(u)  Safety  Considerations:  Because  of  the  use  of  hydrogen  and 

the  desire  to  prevent  premature  firing 
into  low  pressure  accelerator  gas,  the 
system  had  to  have  permissive  and 
positive  firing  control. 

(j)  Energy  Source:  From  existing  VKF  energy  storage  systems 

This  performance  capability  required  a  very  versatile  energy 
storage  system.  After  studying  existing  VKF  energy  storage  systems, 
it  was  concluded  that  a  simulated  transmission  line  made  from  an  exist¬ 
ing  10®- joule,  4-kv  capacitor  bank  would  be  preferable  to  existing 
inductive  energy  storage  systems,  and,  therefore,  theoretical  studies 
were  made  to  determine  the  system  characteristics  and  required  values 
for  various  parameters. 


6.2  THEORETICAL  STUDIES 

Several  case  studies  were  made  to  determine  the  effects  of  various 
parameters.  These  studies  were  made  using  an  analog  computer 
programmed  with  the  exact  differential  equations.  The  following  rela¬ 
tionships  were  used  from  Ref.  39. 

(1)  Characteristic  impedance  for  a  lossless  line 

Zo  =  ohm9 
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The  terms  Ljj  and  Cj^  are  the  inductance  and  capacitance  of  each  sec¬ 
tion  measured  in  henries  and  farads,  respectively. 

(2)  Pulse  Duration 

t  =  2N  [  LnCn  ]‘/j  sec  (24) 

where  N  =  number  of  transmission  line  sections  employed. 

(3)  Load  Voltage 

-  v«  (i  -  -afi-zr)  <25) 

V0  =  initial  voltage  on  capacitors  in  volts 
Rl  =  resistance  of  load  in  ohms 


(4)  Load  Current 

!l  =  vl/rl  amP 


(26) 


In  addition,  the  MHD  accelerator  channel  resistance  was  assumed 
to  be 


(5)  Resistance  of  the  Channel 

R  -  -1  4“  (B  =  0) 

o  At 


(27) 


where  / a  =  effective  distance  between 

accelerator  electrodes,  m 

Ae  =  area  of  accelerator  electrode,  m2 

a  =  electrical  conductivity,  mho/m 


From  the  accelerator  dimensions  given  in  Section  VIII  and  the 
geometrical  considerations  of  Fishman  (Ref.  40),  this  equation  becomes 

R,  =  — 6--x-  - ■■  (B  =  0)  (28) 

£7  La 

where  a  =  electrical  conductivity,  mho/m 

La  =  length  of  accelerator  electrodes,  in. 

(6)  Current  Density  (Neglecting  the  UB  Term) 

j  =  *E  (29) 

where  E  =  electric  field  intensity,  v/m 
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The  following  cases  were  considered: 

6.2.1  Cose  One 

The  simplest  case  considered  was  the  lossless  transmission  line  as 
shown  in  Fig.  24.  For  this  case,  with  a  matched  load  {i.  e. ,  H^/Z0  =  1) 
and  values  of  the  various  parameters  in  the  range  of  interest,  the  effect 
of  varying  the  number  of  sections  in  the  line  {Fig.  25)  was  to  decrease 
the  pulse  duration  as  the  number  of  loops  decreased  while  the  average 
current  level  remained  the  same.  Figures  26  and  27  show  the  behavior 
for  the  matched  case  of  the  loop  currents,  and  the  voltage  across  each 
section  capacitance  as  a  function  of  time.  Figure  28  shows  typical  plots 
of  the  shape  of  the  current  pulse  through  the  load  for  various  values  of 
Rl/Z0  where  the  characteristic  impedance,  ZQ,  was  constant  and  the 
resistance,  R^,  was  varied  according  to  Eq.  (28)  for  a  number  of  elec¬ 
trode  lengths  and  gas  conductivities.  Figures  29,  30,  and  31  summarize 
the  plots  of  Fig.  28  plus  many  others  not  presented  herein  by  showing  the 
variation  of  load  current,  current  density,  and  electric  field  intensity 
with  accelerator  electrode  length  and  gas  conductivity, 

6.2.2  Case  Two 


Having  seen  in  Section  6.  2.  1  that  a  reasonably  square  pulse  could  be 
obtained  across  the  load  with  the  desired  range  of  current  density,  the 
effect  of  adding  inductance  to  the  load  was  considered  as  shown  sche¬ 
matically  in  Fig.  32.  With  the  resistance  of  the  channel,  Rl>  equal  to 
the  characteristic  impedance,  ZQ,  it  was  evident  (Fig.  33a)  that  for  a 
small  value  of  inductance  (10"®  henry)  the  average  current  across  the 
load  remained  the  same  as  for  the  case  with  no  inductance  (Fig.  25)  and 
that  the  inductance  removed  most  of  the  ringing  effect  at  the  initial  part 
of  the  current  pulse.  This  value  of  inductance  changed  the  rise  time  by 
less  than  2  percent.  However,  for  larger  values  of  inductance  (5xl0“® 
and  10“ ^  henries.  Figs.  33b  and  c),  the  rise  time  changed  considerably, 
the  pulse  shape  was  no  longer  square,  and  the  average  current  did  not 
reasonably  approach  the  value  in  Fig.  25.  It  was  therefore  desirable  to 
limit  the  inductance  in  the  load  to  approximately  the  value  of  Ljq  to  obtain 
a  smooth  pulse. 

6.2.3  Case  Three 

The  circuit  in  Fig.  34  considered  the  effect  of  resistance  in  the 
cables  connecting  the  capacitors,  and  calculated  results  for  such  a  cir¬ 
cuit  are  shown  in  Fig.  35.  A  comparison  of  the  pulse  shapes  (Figs.  33a 
and  35a)  revealed  that  a  small  resistance  in  each  loop  (-  5  x  10” 5  ohms) 
had  very  little  effect  on  the  shape  or  value  of  the  current  across  the  load. 
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but  as  this  value  was  increased  to  5  x  10“  3  ohms  (Fig.  35d)  the  initial 
slope  decreased  slightly,  there  was  a  slight  decrease  in  the  average 
current,  and  the  current  began  falling  off  sooner.  It  was  therefore  de¬ 
sirable  to  limit  the  resistance  in  each  section  to  one  milliohm  or  less. 

6.2.4  Cose  Four 


The  final  case  considered  the  total  effect  of  line  and  load  inductance 
and  line  resistance  as  shown  schematically  in  Fig.  36.  It  was  desirable 
to  study  the  current  which  could  be  obtained  for  several  pulse  durations 
as  determined  by  Eq.  (24).  Since  the  value  of  Cj^  could  not  be  decreased 
without  losing  energy  for  a  given  load  voltage  requirement,  five  arbitrary 
values  of  Lj^j  were  selected  to  give  pulse  widths  from  0.  4  to  2.  0  msec  in 
0.4-msec  increments. 

Figure  37  shows  the  behavior  of  the  current  across  the  load  for  one 
of  the  five  values  of  inductance  for  the  matched  case.  This  plot  and  other 
matched  and  unmatched  cases  not  presented  herein  are  summarized  for 
one  of  the  five  values  of  inductance  in  Fig.  38.  The  effect  of  accelerator 
electrode  length  and  gas  conductivity  on  the  load  current  is  given  in 
Fig.  38  and  the  effect  on  current  density  in  Fig.  39. 

Although  all  the  pulse  durations  are  not  presented  herein,  the  pulse 
duration  lengthened  as  the  value  of  inductance  increased  (Eq.  (24)),  but 
as  the  inductance  increased  the  maximum  current  which  could  be  ob¬ 
tained  decreased.  Also,  as  the  inductance  increased,  the  current  density 
which  could  be  obtained  with  a  specific  load  resistance  decreased.  There¬ 
fore,  when  the  current  pulse  lengthened,  the  current  density  obtainable 
decreased. 

Maximum  energy  transfer  in  the  specified  pulse  time  was  obtained 
with  a  matched  (Rj^  =  ZQ)  load.  Figure  40  shows  normalized  energy 
delivered  in  the  first  pulse  as  a  function  of  the  degree  of  mismatch. 
Normalized  current  through  and  voltage  across  the  load  are  also  plotted 
as  a  function  of  normalized  load  resistance  in  Fig.  40. 


6.3  SYSTEM  DESCRIPTION 

Based  on  the  analog  computer  studies  presented  herein  and  other 
local  design  considerations,  the  electric  field  system  configuration  that 
was  built  is  that  which  is  shown  schematically  in  Fig.  41.  Referring  to 
Fig.  41,  the  electric  field  (E-field)  system  consisted  of  a  ten-section, 
simulated  transmission  line,  an  air  gap  switch,  the  load  bus,  and  the 
accelerator  electrodes  through  which  (when  B  =  0)  a  resistive  gas  flowed. 
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Each  of  the  ten  sections  of  the  simulated  transmission  line  con¬ 
sisted  of  capacitance,  line  resistance  and  inductance,  and  a  variable 
inductance.  The  section  capacitance  was  provided  by  50  "can- stacks” 
connected  in  parallel.  Each  can- stack  consisted  of  two  125-«f,  4-kv 
rated  capacitors  connected  in  series  to  provide  62.  5  /^f  at  8  kv  in  each 
can- stack,  making  a  total  section  capacitance  of  3.  125  millifarads  at 
a  design  rating  of  8  kv.  The  line  resistance  and  inductance  was  due 
to  the  parallel  combination  of  200  RG8A/U  coaxial  cables,  each  75  ft 
long,  connecting  the  section  capacitance  to  a  collector  plate.  Each 
can- stack  was  connected  to  the  collector  plate  by  four  parallel  cables. 
The  total  line  resistance  was  1.04  milliohms  and  the  total  line  induct¬ 
ance,  0.0286  microhenry  (ph).  The  75-ft  length  of  each  cable  was  due 
to  a  multiple  accelerator  location  test  plan  and  the  desire  to  keep  the 
inductance,  resistance,  and  capacitance  of  each  section  balanced. 

The  range  of  impedance  matching  as  a  function  of  section  inductance, 
capacitance,  and  total  pulse  time  is  shown  in  Fig.  42,  The  variable 
inductance  was  provided  by  either  1  or  2  turn  inductors  of  variable 
diameter  to  give  the  values  shown  in  Fig.  41.  The  inductors  were 
physically  located  between  the  collector  plates  of  each  section  and 
were  made  from  1-  by  1-in.  cross-section,  heat-treated  aluminum 
bar  stock.  Figure  43  shows  photographs  of  the  racked  capacitors  of 
each  section,  the  current  collector  plate,  and  the  installed  inductors 
between  each  section. 

The  ten  sections  were  connected  in  series  as  shown  in  Fig.  41  and 
were  switched  on  by  means  of  an  air  gap  switch  shown  in  Fig.  44.  The 
current  pulse  was  fed  through  a  coaxial  bus  system*  to  the  accelerator 
electrodes  as  shown  in  the  photograph  in  Fig.  45.  A  discussion  of  the 
accelerator  is  provided  in  Section  VIII. 

Because  of  the  requirement  for  quickly  disconnecting  the  load  bus 
from  the  accelerator  electrodes  and  the  time  element  involved  in  the 
original  fabrication  of  the  load  bus,  a  higher  inductance  (L^us  ^  0.8  Ath) 
was  obtained  than  was  desirable. 

The  "switching  on”  of  the  electric  field  was  accomplished  by  the 
circuitry  shown  schematically  in  Fig.  44.  For  normal  operation,  the 
tungsten-tipped  gap  electrodes  were  spaced  such  that  the  gap  would  hold 
off  2  kv  more  than  the  bank  charge  voltage  to  prevent  premature  self¬ 
firing  caused  by  changes  in  humidity,  dust  in  the  gap,  et  cetera.  The 


*An  improved  bus  system  was  installed  later  because  the  bus 
system  shown  in  the  photograph  was  found  to  transmit  forces  to  the 
electrodes. 
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automatic  gap  firing  system,  designed  by  Mr.  D.  S.  Bynum,  ARO,  Inc., 
sensed  the  pressure  rise  behind  the  shock  wave  with  a  pressure  trans¬ 
ducer  immediately  downstream  of  the  driver  diaphragm  and,  provided, 
the  pressure  level  was  at  a  predetermined  level,  transmitted  a  10-msec 
permissive  pulse,  only  during  which  the  gap  could  be  fired  automatically. 
With  the  permissive  pulse  provided,  when  the  shock  wave  arrived  at  a 
given  station  upstream  of  the  accelerator  electrodes,  the  system  sensed 
the  passing  shock  wave  with  an  ionization  probe  and  transmitted  a  pulse 
to  a  time  delay  unit.  After  a  predetermined  time  delay,  the  third  elec¬ 
trode  in  the  auxiliary  gap  was  pulled  to  ground  potential  via  a  thyratron 
firing  and  the  auxiliary  gap  then  conducted  current.  This  placed  the 
20-kv  auxiliary  capacitors  across  the  gap  switch.  The  very  short  pulse 
from  the  auxiliary  capacitors  broke  down  the  gap,  and  the  voltage 
balance  network  drew  a  small  sustaining  current  until  the  discharge 
across  the  electrodes  was  underway.  Inductance  between  the  gap  and 
capacitors  and  between  the  gap  and  the  accelerator  was  essential  for 
this  mode  of  operation.  The  time  resolution  in  synchronizing  the  switch¬ 
ing  on  the  electric  field  with  the  shock  wave  is  on  the  order  of  a  few 
microseconds.  The  accelerator  electrode  voltage  balancing  circuit 
shown  in  Fig.  46  was  used  to  balance  the  voltage  between  each  electrode 
and  the  accelerator  housing  (ground  potential)  to  relieve  a  critical 
electrode  insulation  problem  inside  the  accelerator. 

The  remaining  components  of  the  electric  field  system  consisted 
of  a  standard  8-kv  charger  and  high  voltage  relays  for  disconnecting 
during  discharge,  a  grounding  switch  system  for  grounding  when  not  in 
use,  a  crowbar  switch,  and  a  crowbar  load  for  dumping  the  energy  in 
the  bank. 

Standard  instrumentation  measurements  included  total  current, 
total  electrode  voltage,  voltage  between  each  electrode  and  the  accel¬ 
erator  housing  (ground  potential),  voltage  across  the  air  gap  switch, 
and  rate  of  change  of  current.  All  voltage  measurements  were  made 
using  a  standard  voltage  divider  and  were  displayed  on  an  oscilloscope 
and  the  oscilloscope  traces  photographed.  The  rate  of  change  of  cur¬ 
rent  was  measured  by  means  of  a  calibrated  search  coil  wound  toroi- 
dally  around  the  inner  conductor  of  the  coaxial  load  bus.  The  current 
measurement  was  made  by  passively  integrating  the  rat e- of -change - 
of-current  measurement.  Both  dl/dt  and  the  current  were  displayed 
on  oscilloscopes  and  the  oscilloscope  traces  photographed. 


6.4  SYSTEM  PERFORMANCE 

To  determine  the  performance  of  the  E-field  system,  about  sixty 
discharges  were  made  into  dummy  load  resistors  which  were  matched 
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to  the  characteristic  impedance  of  the  system.  The  dummy  load  dis¬ 
charges  were  made  at  various  initial  voltage  levels  and  with  different 
inductors.  The  discharge  characteristics  were  in  good  agreement 
with  theory  except  for  the  400- ^sec  case  for  which  the  stray  induct¬ 
ance  of  the  system  was  significant  relative  to  the  inductance  of  the 
inductors.  This  resulted  in  a  longer  pulse  duration  and  a  slower  initial 
rise  time. 

Oscilloscope  traces  of  total  current,  total  electrode  voltage, 
voltage  between  each  electrode  and  ground  (balancing  network),  voltage 
across  air  gap  switch  and  rate  of  change  of  current  for  a  typical  low 
level  dummy  load  run  are  shown  in  Fig.  47. 


SECTION  VII 

MAGNETIC  FIELD  SYSTEM 


7.1  INTRODUCTION 

As  the  MHD  parametric  studies  progressed,  it  became  more  and 
more  evident  that  high-level,  transient,  magnetic  fields  over  large 
volumes  would  be  required  for  high  performance,  pulsed  MHD  accel¬ 
erators.  In  addition  to  the  high  magnetic  field  levels,  the  magnetic 
field  had  to  be  tailored  along  the  channel  axis  at  the  entrance  and  exit 
of  the  accelerator  to  prevent  large  kinetic  energy  losses  in  the  working 
gas. 


As  a  result,  a  study  of  previous,  high  magnetic  field  work  was 
made  with  particular  emphasis  on  the  following  general  magnetic  field 
system  design  criteria: 

1.  Magnetic  Field  Level:  0-10  weber/m^ 

2.  Field  Distribution:  Uniform  throughout  accelerator  electrode 

length  and  decreasing  uniformly  to  zero 
along  channel  axis  in  not  less  than  5  tube 
diameters  (-8  in.  )  at  both  entrance  and 
exit  of  accelerator 

3.  Pulse  Time  Duration:  Minimum  of  2000  f*sec 

4.  Pulse  Shape:  Uniform  through  pulse  time  duration 

5.  Pulse  Synchronization:  Magnetic  field  to  be  on  when  electric 
field  "switched  on"  upon  arrival  of  shock  wave 
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6.  Field  Volume:  A  minimum  consistent  with  the  required  field 

distribution  along  the  channel  axis 

7.  Energy  Source:  From  existing  AEDC-VKF  energy  storage 

systems 

To  meet  the  above  design  criteria,  a  detailed  search  of  the  litera¬ 
ture  was  made  to  determine  the  theoretical  considerations  required  and 
to  compile  the  fundamental  coil  parameter  relationships  necessary  for 
a  coil  design. 


7.2  THEORETICAL  CONSIDERATIONS 

7.2.1  Summary  of  Applicable  Work  by  Others 

In  the  literature  {Refs.  40-68),  coils  for  producing  high  magnetic 
fields  are  generally  classified  into  three  broad  categories:  (1)  d-c 
coils,  (2)  pulse  coils,  and  (3)  superconducting  coils.  The  coils  may  or 
may  not  be  operated  continuously,  and  may  or  may  not  be  cooled,  but 
are  generally  operated  for  longer  durations  of  time  than  pulse  coils  and 
are  usually  cooled  in  some  manner.  Pulse  coils  are  usually  used  in 
applications  where  the  duration  of  the  pulse  time  is  on  the  order  of  milli¬ 
second  or  less  and  almost  always  are  not  cooled  since  during  such  short 
times  the  total  heat  that  can  be  transferred  is  usually  small.  This  pro¬ 
vides  an  inherent  temperature  limitation.  Superconducting  coils  are 
operated  at  temperatures  below  their  superconducting  critical  tempera¬ 
ture  and  can  be  pulsed  or  operated  indefinitely  without  generating 
Joulean  heat.  Because  of  the  availability  of  low  impedance,  unipolar 
generators,  superconducting  and  short  time  pulse  coils  were  eliminated 
very  early  in  the  study.  In  this  section,  the  fundamental  parametric 
relationships  applicable  to  uncooled,  long  pulse  duration  (-  10  sec)  coils 
are  presented. 

7.2. 1.1  Summary  of  Existing  Coils 

Powerful  d-c  coils  have  long  been  used  to  study  various  magnetic 
effects;  however,  by  far  the  majority  of  the  cases  were  such  that  the 
higher  fields  were  produced  over  very  small  volumes.  This  was  un¬ 
fortunate  from  the  MHD  point  of  view  since  high  fields  over  large 
volumes  were  desirable;  however,  much  was  gained  from  studying  the 
work  of  others,  while  at  the  same  time  keeping  in  mind  the  large  volume 
requirements. 

The  first  powerful  d-c  coil  was  constructed  by  Francis  Bitter  at 
MIT  in  1939  (Ref.  69).  The  coil  consisted  of  a  stack  of  flat  copper 
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washers  and  insulating  spacers,  split  and  connected  so  that  the  washers 
form  a  continuous  helix  for  the  current  path.  This  coil  produced 
10  weber /m^  over  a  volume  of25x  10'^m^  and  was  water  cooled. 

For  a  long  time  the  Bitter  coil  was  the  only  one  in  the  world  which 
could  produce  such  high  fields.  In  the  last  decade,  however,  powerful 
coils  have  been  constructed  in  several  countries  and  are  summarized 
in  Table  I  along  with  other  pertinent  information.  In  their  complex  of 
equipment,  cooling  principle  and  design,  they  are  basically  similar  to 
Bitter's  coil. 

The  first  pulse  coil  for  producing  strong  transient  magnetic  fields 
was  constructed  by  Kapitsa  (Ref.  42)  in  1923-1927.  With  this  device, 
magnetic  fields  of  50  weber/m^  were  achieved.  Pulse  techniques  were 
further  developed  in  1956-1961  such  that  regular  investigations  could 
be  made  at  field  levels  of  70  weber /m2.  Since  then,  extremely  high 
magnetic  fields  have  been  generated,  the  highest  reported  in  the  litera¬ 
ture  being  1400  weber/m^,  produced  by  chemical  implosion  at  Los 
Alamos  (Ref.  43).  Table  II  summarizes  the  more  important  of  the  exist¬ 
ing  pulse  coils  along  with  other  pertinent  data. 

From  the  d-c  coils  field-power  relations,  very  high  magnetic  fields 
require  very  high  power  levels.  The  principal  advantage  of  a  pulse 
system  is  the  ability  to  store  energy  over  a  long  period  of  time  and 
discharge  the  energy  through  the  coil  in  a  short  time,  thus  generating 
very  high  transient  power  levels,  and  relying  on  thermal  inertia  to  pre¬ 
vent  overheating. 

Thus,  by  this  crude  definition  of  a  pulse  coil,  the  same  general 
d-c  coil  parameter  relationships  hold  provided  "skin  effect'  is  taken 
into  consideration  when  defining  the  current  density  distribution  in  the 
conductors  of  the  pulse  coil. 

7. 2. 1.2  Coil  Parameter  Relationships 

According  to  Montgomery  (Ref.  70),  of  the  several  ways  of  gener¬ 
ating  a  magnetic  field,  the  choice  is  dictated  by  considerations  of  field 
magnitude,  volume,  duration,  power,  weight,  and  space.  Thus,  it  was 
of  interest  to  summarize  the  relationships  between  field,  current, 
power,  and  volume  for  a  uniform  current  distribution,  a  Bitter  radial 
current  distribution,  a  Gaume  current  distribution,  and  a  Kelvin  current 
distribution.  These  current  distributions  are  sketched  at  the  top  of 
Table  III.  In  general,  the  relations  summarized  in  Table  III  are  essen¬ 
tially  the  same  as  those  of  Montgomery  (Ref.  70). 
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7. 2. 1.3  Magnetic  Stresses  in  Coils 


The  interaction  between  magnetic  fields  and  currents  produces 
forces-  For  a  cylindrical  coil  these  forces  act  in  such  a  manner  as  to 
compress  the  coil  axially  and  enlarge  the  coil  radially.  Thus,  accord¬ 
ing  to  Furth  (Refs.  40  and  41)  the  axial  magnetic  field  of  an  air  core  coil 
behaves  like  a  two-dimensional  gas  that  exerts  a  radial  pressure 
mathematically  identical  with  the  magnetic  energy  density, 

P„  =  newton /mJ 

2 (30) 

=  57.5  B3  lb/in,2 

where  B  =  weber/ma 

=  4it  x  10  7  for  free  space 


Generally,  for  a  thick-walled  coil,  the  stresses  are  simply  treated 
using  thick-walled  pressure  vessel  theory.  Referring  to  the  coil  model 
shown  in  Fig.  48,  Furth  gives  the  radial  stress  at  any  point  (a)  as 


Sa 


where  a 


a2/ai 


(31) 


and  the  tensile  stress  at  any  point  (a)  as 


for  a  large 


(32) 


(33) 

(34) 


Thus,  the  maximum  tensile  stress  is  equal  to  Pm,  Figure  49  shows  the 
approximate  magnetic  yield  points  of  various  metals  in  terms  of  the 
B-field.  Another  way  of  interpreting  Fig.  49  is  in  terms  of  the  highest 
steady- state  magnetic  field  that  can  be  attained  using  various  coil 
materials  before  the  material  begins  to  yield.  Furth  (Ref.  41)  has  made 
a  large  number  of  tests  to  verify  the  above  relations  and  the  attainable 
levels  given  in  Fig.  49. 


The  above  relations  only  relate  to  the  bursting  stresses,  and  do  not 
take  into  account  inertial  effects.  There  are  also  axial  compressive 
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stresses  which  must  be  considered.  According  to  Montgomery  (Ref.  70), 
the  maximum  average  value  of  the  axial  compressive  stress  is  approxi¬ 
mately  20  percent  of  the  bursting  stresses  and  except  for  insulation 
considerations  can  generally  be  ignored. 

7. 2.1. 4  Skin  Effect  Considerations 


In  Bitter-type  pulse  coils,  the  current  distribution  is  governed  by 
skin  depth  phenomena.  The  skin  depth  is  defined  as  the  depth  below  the 
surface  at  which  the  current  density  is  1/e  of  its  value  at  the  surface. 
The  skin  depth,  S,  is  a  function  of  frequency,  resistivity,  and  permea¬ 
bility.  Stratton  (Ref.  44)  gives  the  relation 


where 


P  =  resistivity  in  ohm-meter 

oi  =  frequency  in  radians  per  second  (=  2nv) 

v  =  frequency  in  cycles  per  second 

P  =  permeability  of  material,  4w  x  lO-"7  for 
a  nonmagnetic  medium 


(35a) 


For  copper  at  293°K,  p  =  1.  72  x  10"8  ohm-meter.  Thus,  Eq.  (35a) 
becomes  for  copper 

„  2  *•  1/ 

S  -  (6.6  x  10  }  v  '2  meters  (35b) 

For  a  1-msec  pulse  (-  500  cycle/sec),  the  skin  depth  in  copper  is  only 
0.295  x  10-2  m  or  about  3  mm.  Thus,  it  is  seen  that  most  disc-type 
pulse  coils  have  somewhat  large  unused  portions  of  conducting  medium 
which  are  used  only  for  additional  strength  and  as  a  sink.  If  this  addi¬ 
tional  strength  and  heat  capacity  are  not  needed,  then  it  is  obvious  that 
the  maximum  radial  thickness  of  the  individual  coil  turn  should  be  on 
the  same  order  as  the  skin  depth. 


According  to  Montgomery  (Ref.  70),  in  determining  the  various 
parameter  relations  for  a  pulse  coil,  if  S/a j  >  1,  j  (a)  quite  closely 
approximates  the  d-c  case  of  a  l/a  current  distribution;  however,  in 
most  pulse  coils  0  <  S/ai  <  1  and  the  current  drops  off  much  faster 
than  1/a.  Thus,  if  a  is  close  to  1,  a  very  reasonable  approximation 
can  be  found  by  assuming  the  current  is  l/a  since  the  field  is  not  very 
dependent  on  the  current  distribution  for  this  case.  Further,  it  is 
generally  conventional  to  assume  the  current  is  l/a  when  calculating 
inductance  and  resistance,  and  thus  the  actual  values  of  L  and  L/R 
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can  be  expected  to  be  smaller  and  R  larger  when  assuming  a  1/a  current 
distribution. 

7.2.2  High  Transient  Magnetic  Fields  for  MHD  Accelerators 


For  high  performance  pulsed  MHD  accelerators,  high  transient 
magnetic  fields  are  required  over  large  volumes.  In  addition,  the 
magnetic  field  must  be  tailored  along  the  channel  axis  at  the  entrance 
and  exit  of  the  accelerator  to  prevent  large  kinetic  energy  losses 
(Ref.  19)  in  the  working  gas. 

7.2.2. 1  Type  of  Coil  for  a  Pulsed  MHD  Accelerator 

The  desired  duration  of  the  j  x  B  force  acting  on  the  gas  in  the  high 
performance  pulsed  MHD  accelerator  is  on  the  order  of  one  millisecond. 
Therefore,  the  simplest  type  of  coil  to  use  for  the  B-field  would  be  a 
pulse  coil.  This  is  particularly  true  if  the  switching  synchronization  of 
the  E-  and  B-fields  can  be  readily  achieved.  The  pulse  coil  is  most  de¬ 
sirable  mainly  because  the  complex  cooling  system  can  be  eliminated 
and  thus  a  greater  degree  of  design  freedom  is  allowed  in  the  coil  and 
accelerator  configuration.  However,  skin  effect  presents  a  more  dif¬ 
ficult  situation  in  the  B-field  tailoring  problem.  This  can  be  overcome 
by  providing  for  a  ''forced"  current  density  distribution  in  the  coil  wall 
and  will  be  discussed  more  thoroughly  later. 

In  any  event,  the  superconducting  coil  is  presently  eliminated  if 
the  B-field  levels  required  exceed  about  10  weber/nr.  Although  a  plain, 
cooled  coil  can  be  pulsed  and  has  the  advantage  of  lower  I^R  heat  genera¬ 
tion  since  R  is  lower,  the  overall  gain  is  usually  negative  because  of  the 
refrigeration  system  that  is  reqxured.  Thus,  for  these  reasons,  an 
attempt  should  be  made  to  design  the  B-field  coil  to  operate  in  a  pulse 
duration  regime  where  the  skin  effect  is  low  enough  to  be  tolerated  and  yet 
the  cooling  system  can  be  left  out  of  the  design. 

7. 2. 2. 2  Coil  Configurations  for  MHD  Accelerator 

In  general,  just  about  any  reasonable  coil  configuration  can  be  built 
to  provide  the  B-field  for  an  MHD  accelerator,  provided  funds  are 
available.  The  coil  parameter  relations  presented  in  Table  III  apply; 
however,  an  "effective"  a,  B,  and  a^  will  have  to  be  applied  where  the 
coil  geometry  used  is  other  than  that  described  herein. 

7. 2. 2. 3  Largo  B-Field  Volumes  for  MHD  Accelerator 

Although  in  most  of  the  experimental  work  noted  herein  small 
B-field  volumes  were  used,  large  volumes  are  strictly  a  matter  of 
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funds  (Ref.  51,  i.  e.  ,  coil  and  power  supply  costs).  Since  magnetic 
pressure  is  a  function  of  B^,  a  larger  inside-diameter-coil  presents 
a  greater  problem  in  the  design  of  the  de  facto  pressure  vessel.  A 
plot  of  the  magnetic  pressure  developed,  using  the  relation  (Eq.  (30)) 


2^0 


newton/m* 


or 


Pm  =  57.5  B 1  lb/W 

where  B  is  in  weber/m^,  is  shown  in  Fig.  50  as  a  function  of  B. 


7. 2. 2. 4  Energy  in  the  B- Field  for  MHD  Accelerator 


The  energy  density,  U^,  in  a  magnetic  field  is 

Ud  =  joule/m*  <36) 

To  find  the  total  energy,  one  must  know  B(x,  y,  z)  over  all  space  and 
integrate  over  all  space.  This  is  most  often  difficult  to  do.  By 
thinking  in  terms  of  an  ideal  coil  (one  where  no  flux  linkages  exist 
within  the  coil  wall),  then  every  flux  line  outside  the  coil  must  also 
continue  inside  the  coil.  Therefore,  by  taking  the  inside  volume  (Vj) 
of  the  coil  and  multiplying  by  B0  at  the  center  of  the  coil  yields,  very 
approximately,  one-half  the  total  energy  required. 


Thus,  the  maximum  total  energy,  Urp,  that  can  be  stored  in  the 
magnetic  field  is 


UT 


joules 


(37) 


A  plot  of  the  total  energy  per  unit  inside  volume  (Vj)  of  the  coil  is 
shown  in  Fig.  51.  In  addition  to  the  energy  stored  in  the  magnetic 
field,  the  Joulean  heat  energy  must  be  added  to  arrive  at  the  total 
energy  required. 


7. 2. 2. 5  Attainable  Magnetic  Field  Levels  for  MHD  Accelerator 


From  the  graph  in  Fig.  49,  it  is  seen  that  magnetic  field  levels 
up  to  about  60  weber/m^  can  be  attained  without  destroying  or 
damaging  the  coil.  The  limiting  factor  in  this  case  is  the  tensile 
yield  strength  of  the  coil  material  (assuming  adequate  cooling  or  mass 
to  limit  the  temperature  rise  in  the  coil).  Further,  this  assumes  a 
suitable  structure  can  be  designed  to  contain  the  de  facto  pressure 
vessel.  The  design  of  a  d-c  coil  to  produce  a  steady- state  magnetic 
field  level  of  60  weber/m^  would  be  a  formidable  engineering  problem 
and  would  be  very  expensive. 
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7. 2. 2. 6  Tailoring  the  Magnetic  Field  at  the  Entrance  and  Exit  of  a  MHD  Accelerator 

According  to  Leonard  {Ref.  19)  it  is  necessary  to  tailor  the  mag¬ 
netic  field  at  the  entrance  and  exit  of  a  MHD  accelerator  to  prevent 
high  kinetic  energy  losses  in  the  working  fluid.  This  can  be  done 
simply  by  choosing  the  proper  wall  thickness  for  the  coil  at  the  entrance 
and  exit  of  the  MHD  accelerator.  Figure  52  shows  a  nondimensionalized 
plot  of  the  radial  variation  of  the  axial  field  component  at  several  axial 
stations.  It  should  be  pointed  out  that  this  plot  is  for  constant  current 
density  with  radius,  a  given  a  and  B,  and  it  does  not  consider  skin  effect. 
For  pulse  coils  the  skin  effect  can  be  greatly  reduced  by  a  forced  current 
distribution  through  the  coil  wall.  One  way  to  accomplish  this  is  to  wind 
the  turns  of  the  coil  in  a  spiral  so  that  the  current  in  the  inner  turns  must 
also  flow  in  the  outer  turns.  This  will  not  only  force  the  desired  current 
distribution  in  the  wall,  but  also,  according  to  Furth  (Ref.  40),  will  re¬ 
duce  the  magnetic  pressure  at  the  inside  coil  wall  and  the  temperature 
rise  in  the  skin  volume. 

Based  on  the  high  magnetic  field  theory  and  other  considerations  pre¬ 
sented  herein,  it  was  concluded  that 

1.  Although  the  production  of  B-field  levels  is  a  formidable 
engineering  problem,  levels  of  60  weber/m2  (600,  000  gauss) 
could  be  attained  for  use  with  large  MHD  accelerators. 

2.  To  properly  tailor  the  entrance  and  exit  B-field  distribution  along 
the  accelerator  channel  axis,  a  forced  current  distribution  in  the 
coil  wall  must  be  provided  if  the  skin  depth  is  small  relative  to 
the  coil  wall  thickness. 

3.  The  radial  coil  structure  must  be  designed  such  that  it  can  with¬ 
stand  57.5  B2  lb/in.  2  equivalent  gas  pressure  acting  radially 
on  the  inside  surface  of  the  coil  wall  at  a  given  B-field  level  and 
that  the  axial  coil  structure  must  be  capable  of  withstanding  in 
compression  approximately  20  percent  of  the  radial  pressure. 

4.  Although  coils  have  not  been  built  to  produce  60  weber/m2  over 
volumes  on  the  order  of  a  few  cubic  meters,  such  coils  seem 
technically  if  not  economically  feasible. 

5.  Due  consideration  should  be  given  MHD  pulse  generators  to 
drive  such  a  coil  so  that  the  high  costs  of  power  supplies  can 
be  drastically  reduced. 


7.3  DESIGN  OF  MAGNETIC  FIELD  COIL  FOR  RESEARCH  MHD  ACCELERATOR 

Based  on  the  general  magnetic  field  system  design  criteria  and  the 
theoretical  considerations  presented  herein  and  on  other  studies  not 
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herein  presented,  the  following  specific  design  criteria  were  established 
for  the  present  MHD  accelerator  magnetic  field  coil  design: 

1.  Magnetic  Field  Level:  10  weber/m^  maximum  throughout 

accelerator  electrodes 

2.  Field  Distribution:  Uniform  throughout  accelerator 

electrode  length  of  12  in.  ,  which 
dictated  an  inside  coil  diameter  of  18  in.  ; 
and  uniformly  decreasing  to  zero  in  not 
less  than  8  in. ,  which  dictated  a  forced 
current  distribution  in  the  coil  wall. 

Also,  the  coil  wall  thickness  could  not  be 
less  than  8  in. 

3.  Pulse  Shape:  Damped  RLC  circuit  pulse  shape  with 

seconds  of  relatively  little  change  at  peak 
of  pulse. 

4.  Pulse  Time  Duration:  A  maximum  of  10  sec  because  of  the  time 

constant  of  low  voltage  energy  storage 
system  selected  as  energy  source. 

5.  Pulse  Synchronization:  Because  of  the  long  pulse  time  of  B-field 

pulse.  E-field  to  be  fired  on  peak  of 
pulse. 

6.  Field  Volume:  Inside  diameter  of  coil  to  be  sufficiently 

large  to  provide  the  required  channel  axis 
field  distribution  and  also  to  provide  space 
to  connect  E-field  load  bus  to  accelerator 
electrodes. 

7.  Energy  Source:  Existing  AEDC-VKF  Tunnel  F 

100- Megajoule  Inductive  Energy  Storage 
System.  Maximum  voltage  of  90  v  and 
maximum  current  of  10®  amp,  where  both 
do  not  occur  at  the  same  time  because  of 
decay  in  speed  of  motor- generator  sets. 

8.  Temperature  Rise:  Maximum  AT  ”  125°F 

9.  Coil  Shape:  Circular  design  because  of  high  magnetic 

pressures,  and,  further,  coil  shape  must 
be  compatible  with  accelerator  housing. 
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With  the  above  specific  design  criteria,  the  coil  was  originally 
designed  as  indicated  below. 

From  the  relation  (assuming  uniform  current  distribution) 

Bo  =  fio  J  A  [f  (a,  /3)  ]  weber/m* 

where  the  nondimensional  coil  parameters  a  and  0  are  defined  as  shown 
in  Table  III  and  A  is  the  ratio  of  current  conducting  area  to  total  area. 

For 

B0  =  10  weber/m^ 
a,  =  9  in. 

/io  =  4jt  x  10-?  henry/m 

A  =  0.  9  (assumed  from  other  coil  designs) 

J  =  25,  000  amp/ in.  2  (based  on  allowable  total 

temperature  rise  of  125°F 
in  10  sec) 

the  f  (or,  /3  )  was  computed  to  be 

f(a,0)  =  1.0 

For  the  value  of  f  (a,  )  =  1.  0,  using  the  plot  of  normalized  volume  versus 

j9  with  a  the  running  parameter  as  shown  in  Fig.  53,  the  values  of  a  and  £ 
for  a  minimum  volume  in  the  coil  were  selected.  These  values  were 
taken  as 


a  =  2.555  and  0  =  1.667 

which,  when  actually  calculated,  give  a  value  for  f  (a,  ft )  =  1.  06.  This 
6-percent  margin  was  provided  to  compensate  for  a  reduction  in  field 
strength  at  the  center  of  the  coil  that  was  caused  by  holes  through  the 
coil  walls.  The  holes  were  provided  so  that  the  accelerator  tube  could 
extend  through  the  coil  in  an  area  of  maximum  field  strength. 

From  the  selected  values  of  a  and  ./3,  other  coil  parameters  were 
computed  as  follows: 

Total  Coil  Volume: 

Vt  =  2 jt  at3  [£  (a1  -  1  )  ] 

=  42.2  x  105  in.3 
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Total  Current  per  Turn: 

IT  =  (  VT  *P>  J' 

E  c  o  il 

Allowing  65  v  for  the  coil  after  bus  voltage  drop, 
It  -  255,000  amp 
Number  of  Turns: 


N  = 


A-p  A  J 


N  =  37  turns 


Area  per  Turn: 


An  = 


l 


J  A 

An  =  11.34  in.J 
Outside  Diameter: 

2  a2  *=■  2a  at 


=  46.0  in 


Coil  Length: 


2b  =  2/3 
2  b  =  30.0  in. 
Turn  Thickness: 

An 


Ntot.l  2b 

=  0.378  in. 

—  t»j 

copper  ^totai 

=  0.340  in. 


^insulation  ^lotal 

=  0.038  in. 


copper 
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Magnetic  Pressure  Inside  Coil: 


Pm 


Bo2 


,ldial  2  n0 

=  5750  psi 


A  sketch  of  this  coil  design  is  shown  in  Fig.  54.  A  half- scale  model 
of  this  coil  was  built,  and  measurements  were  made  of  the  magnetic  field 
distribution  along  the  MHD  accelerator  channel  axis  to  determine  the 
effect  of  the  accelerator  holes  through  the  coil  walls.  A  photograph  of 
this  model  is  shown  in  Fig.  55,  and  the  measurements  of  the  magnetic 
field  distribution  along  the  MHD  accelerator  channel  axis  are  shown  in 
Fig.  56. 

The  original  design  of  the  coil  shown  in  Fig.  54  specified  turns  made 
of  copper  sheet,  spirally  wound  to  give  37  turns  and  a  hole  bored  through 
the  coil  walls  for  the  accelerator  tube  to  go  through.  Because  of  a  copper 
shortage  and  long  lead  times  required  for  copper  sheet,  an  alternate  de¬ 
sign  was  chosen  where  copper  magnet  wire  was  substituted  for  the  copper 
sheet.  In  addition,  however,  the  magnet  wire  had  to  be  wound  around  the 
accelerator  tube  hole  to  make  the  inner  spirals,  or  pancakes,  effective. 
This  proved  to  be  difficult  to  do,  and  as  a  result  several  design  changes 
were  made  to  make  fabrication  easier.  A  4-in.  -thick  insulator  was  added 
between  the  coil  halves,  and  smaller  spacers  were  added  between  each 
coil  pancake  to  provide  proper  pancake  spacing  so  that  backup  load  bear¬ 
ing  area  between  turns  would  be  sufficient.  These  necessary  design 
changes  resulted  in  other  changes  in  a,  £,  and  the  number  of  turns  to  get 
the  desired  magnetic  field  level  at  the  center  of  the  coil.  The  final  de¬ 
sign  was  done  in  the  same  manner  as  the  first  design  except  using  new 
values  of  o,  N,  and  A  that  were  compatible  with  the  fabrication  tech¬ 
nique.  A  drawing  of  the  final  coil  that  was  built  is  shown  in  Fig.  57,  and 
a  photograph  of  the  coil  installed  is  shown  in  Fig.  58.  The  final  values 
of  the  various  coil  parameters  are  compared  to  the  original  values  in 
Table  IV. 

7.4  SYSTEM  DESCRIPTION 

An  electrical  schematic  of  magnetic  field  system  is  shown  in  Fig.  59. 
The  system  consists  of  the  coil  (final  design)  previously  described,  the 
Tunnel  F  unipolar  generators,  the  interconnecting  bus,  associated  switch- 
gear,  and  the  necessary  control  system. 

A  pictorial  representation  of  the  Tunnel  F  inductance  energy  storage 
system  is  shown  in  Fig.  60,  and  a  description  of  its  normal  operation  and 
performance  characteristics  is  given  in  Ref.  71. 
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The  switch  SW-4  (Fig.  59)  consists  of  eight,  3000-amp-rated,  dis¬ 
connect  switches  ganged  together  and  actuated  by  an  air  cylinder.  The 
switches  are  mounted  on  20-kv  insulators.  The  bus  connecting  the  coil 
to  SW-4  consists  of  42  paralleled  1000-MCM  cables  approximately 
150  ft  long  and  were  designed  to  allow  a  20-v  drop  in  the  bus  at 
280,  000  amp.  Photographs  of  the  SW-4  installation  and  the  terminating 
cables  are  shown  in  Fig.  61.  Overvoltage  protection  for  the  coil  when 
SW-4  interrupts  the  current  after  a  pulse  is  provided  by  a  resistor  and 
a  voltage  breakdown  disc,  both  connected  in  parallel  with  the  coil 
(Fig.  59).  The  control  system  on  the  generators  only  regulates  the 
generator  terminal  voltage  down  to  approximately  0.  75  v,  and  thus, 
when  SW-4  interrupts,  the  resistor  in  parallel  limits  the  coil  voltage 
to  2000  v  or  less.  The  over-voltage  discs  back  up  the  resistors  and 
are  set  to  break  down  at  a  voltage  slightly  higher  than  2000  v. 

A  typical  coil  charging  cycle  is  shown  in  Fig.  62.  The  long  time 
at  maximum  current  through  the  coil  was  necessary  because  of  the  non¬ 
repeatability  of  the  shock  tube  diaphragm  rupture  and  will  be  reduced 
when  better  repeatability  is  achieved. 


7.5  SYSTEM  PERFORMANCE 

The  measured  magnetic  field  levels  inside  the  coil  are  compared 
with  theory  in  Fig.  63,  and  the  agreement  is  considered  satisfactory. 


SECTION  VIII 
BUFFER  ACCELERATOR 


The  buffer  accelerator  was  located  near  the  entrance  to  the  driven 
tube  (accelerator  goes  through  magnetic  field  coil)  as  is  shown  in  the 
drawing  of  the  research  shock  tunnel  system  in  Fig.  64.  Figure  65 
shows  a  photograph  of  the  installed  accelerator  relative  to  the  B-field 
coil  and  other  shock  tube  components.  A  photograph  showing  the  com¬ 
ponents  of  the  accelerator  is  given  in  Fig.  66.  The  details  of  construc¬ 
tion  and  assembly  can  best  be  seen  in  the  isometric  cutaway  drawing 
shown  in  Fig.  67. 

The  inside  diameter  of  the  accelerator  was  made  40  mm  to  reduce 
the  total  energy  required.  The  accelerator  housing  was  made  of  a  non¬ 
magnetic  303  stainless  steel  and  was  designed  to  contain  a  maximum 
pressure  of  10,  000  psi.  A  Plexiglas®  window  with  three  ports  was  pro¬ 
vided  for  viewing  the  electric  discharge.  The  electrodes  through  the 
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accelerator  housing  and  the  electrode  conducting  surfaces  were  made  of 
BeCu  for  strength  and  conductivity.  Except  for  the  epoxy  resin  insula¬ 
tion  on  part  of  the  backside  of  the  conducting  surface,  all  insulators 
inside  the  accelerator  were  nylon. 

Versatility  was  provided  in  the  design  so  that  any  length  of  electrode 
conducting  surface  from  0.  5  to  12  in.  could  be  used.  Also,  the  current 
path  to  the  conducting  surfaces  could  either  be  through  the  upstream  or 
downstream  electrode  position,  depending  on  the  desired  effect. 

During  E- field  discharge  studies,  some  difficulty  was  experienced 
with  the  epoxy  resin  potting  on  the  conducting  surfaces.  Improved  pot¬ 
ting  techniques  and  an  electrode-to- ground  voltage  balancing  circuit 
(Fig.  46}  proved  to  be  a  satisfactory  solution  to  this  insulation  problem. 


SECTION  IX 
OPERATIONS 

The  operational  testing  of  the  equipment  described  in  this  report  has 
alternated  with  the  shock  tube  operation  described  in  Ref.  1. 

The  initial  tests  were  made  over  the  period  from  February  18  to 
May  21,  1965,  and  consisted  of  discharges  without  the  magnetic  field. 
During  these  tests  the  emphasis  was  on  working  out  instrumentation, 
equipment,  and  timing  problems.  These  tests  indicated  difficulties  with 
the  accelerator,  which  was  redesigned  to  provide  more  electrical  insula¬ 
tion.  These  tests  also  provided  a  base  from  which  to  compare  operation 
with  the  magnetic  field.  Operations  with  combined  electric  and  magnetic 
fields  wrere  undertaken  on  September  21,  1965. 

A  segmented- electrode  accelerator  is  currently  being  designed  for 
the  basic  experiments  on  the  nozzle  accelerator.  Initially,  tests  will  be 
conducted  with  the  accelerator  in  the  same  position  as  the  present  buffer 
accelerator,  with  a  complete  nozzle  installation  at  a  later  date  if  the 
initial  results  are  encouraging.  During  the  time  period  between  the  con¬ 
clusion  of  the  buffer  accelerator  tests  and  the  segmented-electrode  tests, 
a  number  of  shock  tube  tests  will  be  made  (Ref.  4). 
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A  E  DC-T  R -66-25 


Test  Section  Velocity  -  30,000  ft/sec 
Density  Altitude  ■  200,000  ft 

Turbulent  Boundary  Layer 
Nozzle  Half-Angle  *  5  deg 
6  -  26* 

(Based  upon  empirical  results  obtained  for 
lower  velocities) 


Fig.  2  Estimate  of  Usable  Core  Section  Diameter 
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Percent  Frozen  Oxygen 
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AEDC-T  R-66-25 


1962  U.  S.  Standard  Atmosphere  (Ref.  13) 
Air  Data:  Humphrey  and  Nee! 


Fig-  4  Altitude  for  Density  and  Frozen  Speed  of  Sound  Matching 
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Fig.  6  Illustration  of  Duplication  Lines 
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Fig.  8  Schematic  of  Accelerator  Processes  in  x-1  Plane 
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Fig.  10  Processes  for  Reduced  Pressure  at  Accelerator  Exit 
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Fig.  12  Shack  Speed  Ratio  across  MHO  Buffer  Accelerator 


13  Schematic  and  Wave  Diagram  for  MHD  Buffer  Accelerator  Shock  Tunnel 
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Pi  atm  (for  ideal  gas  shock 
tube  with  no  attenuation) 
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Fig.  16  Summary  of  Theoretical  Shock  Tube  Performance  for  Various  Energy  Additions 
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Fig.  17  Ratio  of  Stagnation  Enthalpy  to  Static  Enthalpy  behind  a  Moving  Shock  in  Air 
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Fig.  18  Duplication  Capabilities  of  MHD  Buffer-Augmented  Shock  Tunnels 
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Fig.  21  Nozzle  Accelerator  Schematic 
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Fig.  24  Schematic  Diagram  far  a  Lossless  Transmission  Line  —  Case  One 
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Fig.  28  Concluded 
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Fig-  29  Load  Current  as  a  Function  of  Accelerator  Electrode  Length  for  Different  Gas  Conduct!  viti 


Fig.  31  Electric  Field  Intensity  os  a  Function  of  Accelerator  Electrode  Length 
for  Different  Gas  Conductivities 


Fig-  32  Schematic  Diagram  far  Lossless  Transmission  Li 


ith  a  Load  Consisting  of  Resistance  and  Inductance  -  Cose  Two 


AEDC-T  R -64-25 


CD 


O 


=  8000  volts 
o 


(1)  5.0  x  10“  ^ 

(2)  0.5  x  10-3 

(3)  1.0  x  10-"3 

(4)  5.0  x  10“3 


Fig.  34  Schematic  Diagram  for  a  Transmission  Lino  with  Uniform  Resistance  in  Each  Loop  —  Case  Three 
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Fig.  36  Schematic  Diagram  lor  a  Transmission  Line  with  Inductance  ond  Resistance  in  Series  with  the 
Capacitor  of  Each  Loop  —  Case  Four 

CO 

-5 


AEDC-TR-66-25 


Il>  amp  x  10 


A  E  DC  -T  R-66-25 


Fig.  38  Load  Current  as  a  Function  of  Accelerator  Electrode  Length  and  Gas  Conductivity  for  L 
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Subscripts: 

M  -  denotes  matched  case 
L  -  denotes  load 

o  -  denotes  condition  initial  capacitor 


Fig.  40  Nondimensional  Power,  Energy,  Current,  and  Voltage  as  o 
Function  of  Impedance  Matching 
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Fig.  41  Schematic  Diagram  for  Simulated  Transmission  Line  Electric  Field  System  —  Final  Configuration 
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Fig.  42  Range  of  Impedance  Match  as  a  Function  of  ,  and  Total  Test  Time 


103 


104 


Racked  Capacitors  in  One 
Section 


Fig-  43  Photographs  of  E 


E-Field  Simulated  Transmission  Line 
0-8  kv 
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Schematic 


Fig.  44  Schematic  and  Photograph  of  Air  Gap  Switch 
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a.  Schematic 


Electrode 


Without  Balancing  Network 


b.  Typical  Electrode  to  Ground  Voltage  Traces 
Fig-  4 6  Accelerator  Electrode  Voltage  Balancing  Network 
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Fig.  47  Typical  Electrical  Discharge  Characteristics  of  Electric  Field  System 
Using  Dummy  Load 
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#-  Denote  Magnetic  Field  Lines 
Perpendicular  to  Paper 


Fig-  48  Magnetic  Stresses  in  a  Helical  or  Single-Turn  Coil  Caused  by 
Magnetic  Pressure 
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(Based  on  Tensile 
Strength  at  293°K) 


Fig.  49  Maximum  Magnetic  Field  -  Limited  by  Yield  Stress  of 
Coil  Material  -  All  Coils 
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Fig.  50  Magnetic  Pressure  Developed  inside  o  Coil  os  a  Function 
of  Magnetic  Field  Level 
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Fig.  51  Total  Energy  per  Unit  |n$ide-Coil~Yolume  Stored  in  a  Magnetic  Field 
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Fig.  52  Dimensionless  Axial  Field  of  a  Thick  Coil,  a  =  3,  /3  -  2,  J  Constant 
with  Radius 
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Fig.  53  Normalixed  Coil  Volume  ai  o  Function  of  ft  ond  a 
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Fig.  54  Drawing 


Section  A -A 


gn 


First  Coi I  Desi 


118 


□  .  Along  Tube  Axis 

Fig-  58  Photographs  of  Magnetic  Field  Coil 
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b.  Along  Coil  Axis 
Fig.  58  Concluded 
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Fig.  59  Electrical  Schematic  of  Magnetic  Field  System 
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Electrical  equipment  for  this  system 
0  INDUCTION  MOTOR 
0  EDDY-CURRENT  COUPLING 
0  FLYWHEEL 
0  ACYCLIC  GENERATORS 
0  TRANSFER  SWITCHES 
0  STORAGE  COIL 
0  LOAD  ISOLATING  SWITCH 


includes: 

0  SPECIAL  PROTECTIVE  GAP 
0  MOTOR  CONTROL  CENTER 
0  METAL-CLAD  SWITCHGEAR 
0  AMPLIDYNE  MG  SETS 
0  LOAD  CENTER  SUBSTATION 
0  CLOSED-CIRCUIT  TELEVISION 


Fig.  60  Pictorial  Representation  of  Tunnel  F  100-megojoule  Inductive  Energy  Storage  System 
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Accelerator 

Fig.  62  Typical  B-Field  Coil  Charging  -  Discharging  Cycle 
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Fig.  63  Comparison  of  Measured  Magnetic  Field  Distribution  to  Theory  —  Final  Coil  Installation 
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Fig.  63  Concluded 
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Fig.  66 


Photograph  of  Components  of  Buffer  Accelerator 
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Fig  67  Cutaway  Isometric  Drowing  of  Buffer  Accelerator 
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TABLE  t 

SUMMARY  OF  D-C  COILS  IN  OPERATION 


Name  of  Laboratory 

Field 

Strength  i 
weber/m* 

Working 

Opening, 

cm 

Power 

Consumption, 

mw 

Cooling 

Method 

Supply 

Year  of 
Construction 

Principal 

Investigator 

Lit. 

Ref. 

No. 

Unlv.  in  Leiden 

6.0 

8 

1.4 

Orthodichloro- 

benzene 

Rectifier 

1960 

D.  de  Klerk 

52 

Los  Alamos,  USA 

8.0 

6.3 

Liquid  Hydrogen 

1959 

H.  L  Laquer 

53 

Lyons,  France 

7.0 

1.3 

Water 

Motor-Generator, 

7  ka,  190  v 

1958 

54 

Bell,  Murray-Hill, 
USA 

8.8 

5 

1.5 

Water 

Motor-Generator 

1957 

1.  F.  Kunzler 

55 

Univ.  in  Berkeley, 
USA 

n 

10 

7.5 

Kerosene 

Motor-Generator, 
8.5  ka,  700  v 

1959 

W.  F.  Giaque 

56 

Mass.  Inst,  of 

Tech.,  USA 

10.0 

2.5 

1.7 

Water 

Mot  or -Generator, 

10  ka,  170  v 

1939 

F.  Bitter 

69 

Research  Center  in 
Cleveland,  USA 

10.0 

30 

— 

— 

1961 

mw 
£  i 

57 

Unii/-  of  Tokyo 

n 

— 

3 

— 

- 

Under 

Construction 

S.  Maeda 

58 

Nat'l.  Magnetic  Lab., 
Boston,  USA 

12.6 

2.5 

1.88 

Water  in  Local- 
Boiling 

Motor-Generator, 

10  ka,  170  v 

1961 

H.  Kolm 

59 

Nat'l.  Magnetic  Lab. , 
Boston,  USA 

Water  in  Local- 
Boiling 

Four  2-mw 
Motor-Generators 

Under 

Construction 

F.  Bitter 

60 

NASA,  Cleveland 

Ohio 

30 

Liquid  Neon 

Unipolar  Gen. 

76, 000  amp 
at30  v 

Under 

Construction 

J.  C.  Lawrence 

57 
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TABLE  II 

SUMMARY  OF  PULSE  COILS  IN  OPERATION 


Marne  of  Laboratory 

Magnetic 

Field, 

weber/irr 

Working 

Opening, 

cm 

Pulse 

Length, 

sec 

Cooling 

Supply 

Year  of 
Construction 

Principal 

Investigator 

Lit. 

Ref. 

No. 

Phillips,  Holland 

10.0 

- 

2. 5  x  10'6 

- 

Capacitor,  4  kv, 

2  pf 

1961 

W.  DeBlois 

61 

GERN 

20.0 

5 

IQ'2 

Oil 

Capacitor  Bank. 

3x  Ur  joule 

1960 

M.  Morquo 

62 

Max  Planck  Inst. , 
Munich,  Germany 

30.0 

10-5 

Inside  Liquid  N 

Capacitor  Bank, 

3  x  Ur  joule 

1960 

H.  Bergmann 

66 

Harvard  Univ. ,  USA 

60.0 

0.45 

5  x  10'5 

— 

Capacitor  Bank, 

3  kv.  1725  pt 

1956 

H.  Furth 

40 

Moscow  State  Univ. 

65.0 

0.45 

10‘4 

- 

Capacitor  Bank, 

3  kv,  1200  pt 

1957 

V.  R.  Karsik 

63 

Physics  tnst.,  Acad, 
of  Sciences, 

Moscow,  USSR 

70.0 

0.45 

10-4 

Capacitor  Bank, 

5  kv,  1000  pf 

1961 

V.  R.  Karsik 

63 

Mass.  Inst,  ol  Tech, 
USA 

75. 0 

0.45 

io-4 

- 

Capacitor  Bank, 

2  kv,  3000  p  f 

1957 

S.  Foner 

H.  Kolm 

64 

Harvard  Univ.,  USA 

160.0 

0.4 

10*5 

— 

Capacitor  Bank, 

4  kv,  3000 pf 

1957 

H.  Furth 

40 

Los  Alamos,  USA 

1400.0 

0.3 

2  x  10'6 

-- 

Energy  of  Explosion 
and  Capacitor  Bank, 
20  kv.  7.5  pf 

1960 

C.  Fowler 

43 

KOL  Leiden, 

The  Netherlands 

30.0 

0.6 

IQ’2 

Liquid  H2 

Capacitor  Bank, 

5000  joule 

1961 

J.  C.  A. 
van  der  Sluijs 

65 

Atm.  Wpn.  Res.  Est. 
Aldermaston,  Erg. 

10.0 

15.25 

1 

Liquid  Nj 

Bat.  6700  amp  at 

900  v 

1961 

S.  Skellett 

67 
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TABLE  III 

SUMMARY  OF  CURRENT,  FIELD,  AND  POWER  RELATIONS  FOR  DIFFERENT 

CURRENT  DISTRIBUTIONS 


Unllortn  Current  Bitter  Radial 


Distribution  Current  Distribution 


Current  Distribution 


j(a,xl-Jj 


(n  0  +  la?  *  p2*i;? 


10 


+  <g2-rP  ¥'2 
1  +  II  +  p2)1,z 


Field  Factor.  Fta.pi 

L  “ 

FiekI vs  Power  B0"Gto.pi^^) 

Field  vs  Current  Density  B0  ■  Fla,pi[iajX]  x  10”® 


x  10'1 


Jlal 


I8itx  10‘Ztp" 

u? 

in 

a  +  ta2  ♦  P2)1,z 

4tix  10‘2" 

’"fl  P  ♦  <1  *  P2)w 

la2 -11 

_  l  +  tl  +  pi1'2 

P  in  o 

Tp  +  (a2+p2H 

[4ttP  in  a)1,z  Gta.pi] 


B0-Cfa 


H%) 

B0-F(a,pifja1X]  xlO-6 


Field  vs  Ampere- 
Turns 


Resistance 


Bo-Fta-p{^Trjx10'6 

.  ,_N^  irtotl) 

Xaj  (a -lt$ 


Ffa.P) 


~  Ml 

Zpla  -  It 


R .  71 

X3j  p  in  a 


x  10*^ 
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Current  Distribution 

Geometry  Factor,  Gin,  p ) 

Field  Factor,  Fta.p) 

Field  vs  Ffraer 

Field  vs  Current 
Density 

Field  vs  Ampere- 
Turns 

Resistance 


TABLE  Hi  (Concluded) 


Gaume  Current 
Distribution 


Kelvin  Current 
Distribution 


i<a,  *>  *  ii 


a/a. 


(a/aji2  +  y2 


2  13/2 


4n«  10~2 
in  a 


1/2 


-2FI#,K) 


1/2 


tan 


1/2 


Fie,  Kl  ■  Elliptic  ol  1st  Kind 
e  ■  tan'1  P,  K  -  a 2  -  LA* 

[».»'! !! 


B0  •  Gla 


•«(*)" 


B0  •F|o,pi[j  ajAjx  W6 


1/2 


B0-  Fln.pl  J(NII/[zp(o  -  ll]  j  x  10** 


R  ■  la,pM^AI  J  [llOnl/ln  -  ll][G<a,p]  J2 

f,  •  thickness  ol  cenler  lurn  in  meters 
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TABLE  IV 
COIL  DATA 


First  Design 

Final  Design 

Bn  ,  weber/m2 

10.6 

9.91 

°max 

'max'  amf) 

255, 000 

280, 000 

vmax»  v 

66.9 

74.2 

Rmax,  microlim 

262.5 

265 

L,  microhenry 

555 

520 

^max'  amp/in.2 

25, 000 

25, 000 

Number  of  Turns 

37 

39 

Aturn,  in-2 

11.34 

11.2 

ATmax,  °F/10-sec  Pulse 

125 

125 

Coil  Dimensions 

I.D. ,  in. 

18.0 

18.375 

0.  D.,  in.  (Less  Banding) 

46.0 

50.5 

Width,  in. 

30.0 

37.6 

0  3  32^t3j 

2.55 

2. 76 

p  =  b/a j 

1.67 

2.05 

f<M) 

1.06 

1. 086 

X 

0.9 

0.803 

Conductor  Material 

Copper 

Copper 

Conductor  Shape 

Sheet 

Magnet  Wire 

Banding  Material 

Glass  Tape 

Glass  Tape 
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